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Abstract  of  Dissertation  Presented  to  the  Graduate  Council  of 
the  University  of  Florida  in  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of  Doctor  of  Philosophy 
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by 

PATRICK  CARL  KANGAS 
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Chairman:     Howard  T.  Odum 

Major  Department:     Environmental  Engineering  Sciences 

Energy  models  and  evaluations  were  used  to  understand  form  and 
process  in  landscape  systems.     Embodied  energy  was  used  as  a  compara- 
tive measure  of  value  for  natural   land  forms,  ecosystems,  and  human 
influences  on  landscapes. 

Actual  potential  energy,   embodied  energy,  and  energy  transforma- 
tion ratios  were  calculated   for  ten  landforms  using  energy  analysis 
techniques  and  data  from  the  literature.     Included  were  floodplain, 
salt  marsh,  coral  and  oyster  reefs,  desert  sand  dune  and  arroyo ,  marine 
mud  mound,   ice-cored  palsa  mound,  mountain  glacier  and  spoil  mounds 
formed  by  phosphate  strip  mining  in  central  Florida.     Energy  transfor- 
mation ratios  were  high,  ranging  from  10^  solar  Cal/Cal   for  a  spoil 
mound  to  lO1^  solar  Cal/Cal   for  a  coral  reef. 

Embodied  energies  of  ecosystems  overlying  landforms  were  compared 
with  those  of  landforms.     Embodied  energy  is  109  to  1015  solar 
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Cal/m    for  organization  resulting  from  ecological  succession  and 
10^  to  10^  solar  Cal/m^  for  biomass.     The  difference  between  the 
embodied  energy  in  biomass  and  that  in  successional  organization  is 
suggested  as  a  measure  of  embodied  energy  in  ecological  complexity. 
There  is  up  to  1000  times  more  energy  embodied  in  this   invisible  struc- 
ture than  in  biomass  of  the  three  ecosystems  for  which  data  are  avail- 
able. 

Phosphate  mining  and  reclamation  were  studied  with  ecological 
field  measurements  of  succession,  a  simulation  model  of  succession  and 
energy  evaluations.     Data  on  succession  were  presented  for  a  series  of 
old  spoil  mounds  ranging  from  1  to  50  years  postmining.  Measurements 
included  vegetation  cover,  optical  density,   litter  weight  and 
thickness,   soil  organic  matter,  number  of  ant  nests,   and  tree  height, 
density,  basal  area,  and  species  composition.     Patterns  of  development 
of  the  field  parameters  over  time  demonstrated  that  succession  restored 
the  mined  land  in  a  relatively  short  period. 

For  old  mined  lands  that  have  not  been  reclaimed,  a  threshold  age 
of  25  years  was  calculated  beyond  which  more  energy  accumulated  through 
succession  than  mechanical  reclamation  provided.     Leveling  of  mounds  in 
reclamation  may  divert  a  significant  amount  of  energy  rather  than 
restore  value. 

A  simulation  model  of  succession  was  constructed  and  validated 
with  litter  measurements  on  the  spoil  mounds.     In  these  simulations, 
succession  was  accelerated  through  the  addition  of  seeds  and  litter  so 
that  a  mature  forest  developed  in  33  years  rather  than  50  years. 
Recommendations  were  given  for  managing  succession  as  a  more  economical 
alternative  to  conventional  reclamation. 
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CHAPTER  1 
INTRODUCTION 

Stretching  across  the  continents  of  the  earth  is  a  seemingly  endless 
variety  of  land  types.     Each  individual  area  of  land  surface  has  a  geomor- 
phic  form  and  an  ecological  system  symbiotically  coupled.     Although  some- 
times separated  for  study,   these  features  operate  together  in  natural 
economies  driven  by  the  sun  and  other  energy  sources.     The  organization  of 
these  units  is  landscape,  which  is  expressed  as  particular  spatial 
arrangements  or  patterns.     This  is  a  study  of  the  energy  basis  of  various 
aspects  of  the  landscape.     Considered  are  energy  theories  of  landscape 
along  with  classification  and  calculations  of  component  land  types. 
Models  and  energy  analysis  calculations  are  used  to  consider  how  landform, 
ecosystems,  and  land  use  by  humanity  may  interact  as  a  single  system  and 
to  evaluate  the  amount  of  energy  embodied  in  these  systems. 

A  landscape  model  showing  main  energy  storages  and  pathways  is  given 
in  Figure  1.     Using  land  as  a  platform,   local  energies  such  as  sun,  wind, 
rain,  uplift,  and  others  create  form  and  ecosystems  on  any  given  area. 
These  are  long-term  storages  of  landscape.     Because  these  storages  gener- 
ate  feedbacks   to  secure  and  amplify  energy  inputs,   they  have  value.  As 
noted  by  Odum  (1975:147),  geological  structures  of  landforms  "served  as 
quality  upgraded  potential  energy,   information,  and  material  storages  that 
feed  back  control  actions  on  the  energy  input  processes." 
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At  one  scale,  landform  is  a  consumer,  causing  a  convergence  of  ener- 
gies collected  by  surrounding  land  area.     In  return,   the  landform  performs 
services  such  as  organizing  erosion  and  deposition,  vegetation  patterns, 
drainage,  etc.     Many  individual  land  units  such  as  shown  in  Figure  1  are 
combined  into  landscape  by  exchange  with  the  surrounding  systems  mainly 
through  flows  of  water,  wind,  seed,  and  animal  movements.     These  flows 
integrate  and  organize  the  larger  surface,  generating  properties  of  hier- 
archy in  spatial  dimension  and  feedback  mechanisms  that  support  self-main- 
taining stability. 

Some  of  the  actual  stored  energy  in  landform  is  potential  against 
gravity,  either  to  fill  in  for  depression  forms  or  to  erode  away  for  ele- 
vated forms.     Other  actual  energy  is  stored  in  chemical  concentrations  of 
rocks,  minerals,   and  sediments   (Gilliland  et  al.   1978).     Embodied  energy 
in  landform  is  the  energy  that  was  used  to  develop  the  form  even  though 
most  of  the  actual  energy  used  is  no  longer  stored  in  the  structure.  Cal- 
culations of  embodied  energy  stored  in  individual  landform  types  were  made 
by  multiplying  the  total  rate  of  energy  used  by  the  system  by  the  equiva- 
lent basis  of  solar  calories  and  by  the  length  of  time  needed  to  develop 
internal  storages. 

Embodied  energy  of  landforms  may  be  related  to  effects  of  the  form  on 
the  system.     Can  the  amount  of  control  exerted  by  the  structure  be  in- 
ferred from  its  embodied  energy?     Control   is  shown  in  Figure  1  by  the 
feedback  from  landform  to  the  production  process  of  the  surrounding  land- 
scape.    Calculations  were  made  to  determine  which  is  more  controlling,  the 
living  ecosystem  or  the  nonliving  landform. 

In  landscape  dominated  by  humans,   feedback  from  landform  goes  to 
farms,  housing,  and  other  land  uses.     What   is  the  magnitude  of  control  by 
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land  forms  of  natural  and  human  landscapes?     Does  the  embodied  energy  in 
landforms  influence  the  locations  of  developments  such  as  housing,  utili- 
ties, waste  disposal,  etc.?     Are  there  ways  of  making  better  use  of  land- 
form  energy  for  human  economies? 

In  surface  mining,   feedback  from  landforms  is  interrupted,  and  entire 
landscapes  can  be  totally  disrupted,   leaving  new  landforms  and  surfaces  to 
be  repaired  by  successional  work  of  ecosystems  and  reclamation  efforts  by 
humans.     Energy  analysis  was  used  to  study  phosphate  mining  and  changes  in 
landforms  that  result.     What  mining  alternatives  generate  more  useful 
landscapes?     What  reclamation  designs  best  match  the  available  energy 
inputs?     How  much  energy  is  embodied  in  mining  topography,  and  how  does  it 
compare  with  energy  used  for  reclamation?     Can  succession  be  managed  as  a 
viable  reclamation  alternative?    How  do  these  conclusions  affect  recommen- 
dations for  mining  laws  and  techniques  of  reclamation?     Models  and  energy 
measurements  of  phosphate  mining  in  Florida  were  made  to  compare  natural 
and  human  ways  of  managing  landscape. 

Previous  Study 

Energy  and  Landscape 

The  ancient  Chinese  art  of  geomantics  was  a  system  of  landscape 
design  principles  and  techniques  used  for  siting  developments.  Partly 
scientific  and  partly  magical,   the  goal  of  geomantics  was  "to  produce  a 
landscape  which  had  to  preserve  certain  spiritual  values  and  also  to  ful- 
fill the  practical  purpose  of  supporting  a  dense  population"  (Abercrombie 
1947:30).     The  approach  was  to  design  with  energy-  flows  of  earth  spirits. 
Geomantics,   known  as  Feng-Shui  or  literally  "wind  and  water,"  attempted  to 
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read  and  influence  the  lines  of  flow  of  these  energies  by  concentrating  or 
dispersing  them.     Michell  (1975:14)  states,  "as  traditional  Chinese  medi- 
cine by  acupuncture  treats  the  human  body  by  regulating  the  vital  energy 
that  flows  through  the  skin,   so  the  geomancer  treats  the  body  of  the 
earth.     Textbooks  of  Feng-Shui  describe  the  nature  of  the  telluric  cur- 
rent, give  examples  of  the  energy  patterns  it  forms  in  relation  to  differ- 
ent kinds  of  landscape  and  landscape  features,  and  show  ways  of  adapting 
its  natural   flow  to  human  convenience." 

More  recent  applications  of  energetics  to  landscape  have  been  made  in 
several  contexts.     Price  (1947,  1955,   1959)  was  one  of  the  first  to  recog- 
nize the  importance  of  energy  in  the  development  of  landforms.     He  devel- 
oped a  shoreline  classification  with  five  energy  types:  wave,  current, 
river,  tide,  and  sunlight  (where  organic  systems  dominate).    Within  the 
wave,  current,  and  tide  classes  Price  recognized  at  least  two  intensities, 
"zero"  to  low  energy  and  medium  to  high  energy.     Tanner  (1960a,   1960b)  and 
Parker  (1975)  provide  qualitative  maps  relating  intensities  of  coastal 
energies  to  the  resultant  systems.     Coastal  ecological  systems  were  clas- 
sified according  to  characteristic  energy  sources  by  Odum  et  al .  (1974). 
This  classification  included  nearly  50  different  kinds  of  ecosystems  with 
special  combinations  of  energies  that  generate  them.     The  conclusion  to  be 
reached  from  these  studies  is  that  specific  energy  regimes  cause  or  result 
in  specific  landform  and  ecosystem  types. 

The  related  concept  of  entropy  has  been  used  to  explain  the  develop- 
ment of  landscape  properties  such  as  longitudinal   stream  profiles  and 
drainage  network  patterns.     Leopold  and  Langbein  (1962)  explained  the  ob- 
servation that  rivers  tend  to  have  an  equal  distribution  of  energy  degra- 
dation along  their  longitudinal  profile  by  applying  Prigogine's  (1955) 
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concept  of  systems  trending  toward  a  minimum  rate  of  entropy  generation. 
Since  entropy  generation  is  minimized  they  suggest  that  the  work  done  by 
river  systems  also  tends  toward  a  minimum.     Leopold  and  Langbein  call  this 
the  "principle  of  least  work."     In  a  very  different  context,  but  still  at 
the  landscape  level  of  organization,  Zipf  (1949)  independently  has  devel- 
oped a  principle  of  least  effort.     This  deals  with  distance  minimization 
of  human  movements  over  highway  networks.     This  principle  is  contrasted 
with  the  maximum  power  (or  work)   principle  by  Odum  (1983)  who  notes  that 
the  minimum  entropy  concept  only  applies  to  low  energy  situations  and  that 
as  available  energy  increases  systems  that  survive  tend  to  maximize  rather 
than  minimize  useful  work.     In  fact,  the  principles  of  least  work  do  not 
seem  to  apply  to  the  whole  production  process,  rather  only  to  the  antece- 
dent steps.    Especially  in  terms  of  human  movements,  least  effort  as  dis- 
tance minimization  is  an  important  part  of  maximizing  power,   since  it 
deals  with  bringing  the  factors  of  production  together  to  initiate  the 
production  process.     Jaynes  (1979)   proposed  that   landscape  develops  maxi- 
mum entropy.     Yang  (1971a,  1976;  Yang  and  Song  1979a)  makes  extensive  use 
of  a  theory  of  minimum  rate  of  energy  dissipation  derived  from  entropy 
concepts,  similar  to  that  of  Leopold  and  Langbein  (1962),  to  explain 
fluvial  hydraulics  and  stream  morphology.     According  to  this  theory  "an 
alluvial  channel  is  in  an  equilibrium  condition  if  its  rate  of  energy  dis- 
sipation is  at  its  minimum  under  the  given  climatic,  hydrologic, 
hydraulic,  geologic,  and  man-made  constraints.     The  minimum  value  depends 
on  the  constraints  applied  to  the  system.     If  an  alluvial  channel  deviates 
from  its  minimum  value,   it  will  adjust   its  velocity,  slope,  roughness, 
channel  geometry,   and  pattern  in  such  a  manner  that  the  rate  of  energy 
dissipation  can  be  minimized  to  regain  the  equilibrium"  (Yang  and  Song 


1979b:56).    Many  derivations  are  given  by  Yang  illustrating  the  theory 
including  those  showing  that  meanders  (1971b)  and  pool  and  riffle  sequen- 
ces (1971c)   are  means  for  rivers  to  minimize  the  rate  of  energy  dissipa- 
tion.    However,  at  least  with  meanders,   the  river  may  be  subsidizing  the 
floodplain  forest  ecosystem  surrounding  the  channel  to  maximize  energy 
flow  by  developing  rich  bottomland  soils  and  spreading  out  water  to 
increase  its  interaction  with  other  energy  flows.     Forests  growing  under 
these  conditions  often  have  high  gross  primary  productivity  (Brown  et  al . 
1979),  which  indicates  maximum  work.     Apparently  the  question  of  when 
landscapes  maximize  and  minimize  flow  requires  more  study.     In  another 
thermodynamic  analogy  Scheidegger  (1970)  substitutes  elevation  of  land 
above  baseline  for  temperature  and  land  mass  for  heat  storage,   for  use  in 
gas  laws  and  diffusion  equations.     These  analogies  allow  derivation  of 
characteristics  of  slopes  and  mountain  ranges. 

In  other  energy  applications  in  geomorphology ,  Wolman  and  Miller 
(1960)  suggest  that   for  certain  geomorphic   forces,   the  greatest  amount  of 
work  is  accomplished  by  moderate  events  rather  than  infrequent  catastro- 
phies  or  frequent  but  insignificant  events.     Wolman  and  Gerson  (1978)  mak 
a  distinction  between  sediment  transport  work  (quantity  of  material  trans 
ported  over  a  given  distance  in  a  given  period  of  time)   and  landforming 
work  or  the  actual  sculpture  of  the  landscape.     They  go  on  to  discuss  con 
straints  on  the  geomorphic  effectiveness  of  climatic  events  of  various 
magnitudes  and  time  intervals.     Further  support  is  given  for  the  earlier 
statement  that  different  kinds  and  amounts  of  work  are  done  by  different 
magnitudes  and  frequencies  of  forces.     Carson  and  Kirkby  (1972:72)  provid 
a  table  of  the  orders  of  magnitude  of  geomporphic   forces  and  their  effici 
ency  in  transporting  debris.     Included  are  estimates  for  gravity,  water 
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flow,  rainfall  impact,  and  expansion  forces  due  to  changes  in  moisture, 
freezing,  and  temperature  changes.     They  note  "although  very  large  amounts 
of  thermal  and  biochemical  energy  are  applied  to  the  landscape,  only 
minute  amounts  of  these  forces  are  actually  effective   in  transporting 
debris."     Water  flow  is  suggested  as  the  most   important  geomorphic  force 
both  in  terms  of  energy  expended  and  in  amount  of  material  transport. 
Systems  diagrams  of  streams  are  given  by  Bull  (1979)  showing  feedback 
relationships  between  variables  that  result  in  oscillations  of  sedimenta- 
tion and  erosion.     Streams  are  viewed  as  sediment-transporting  machines 
and  are  analyzed  in  terms  of  the  availability  of  stream  power  to  do  work. 
This  analysis  centers  around  the  concept  of  critical  power  threshold, 
which  is  the  ratio  of  available  stream  power  to  the  critical  power  needed 
to  erode  and  transport  sediments.     Other  geomorphology  texts  also  include 
introductory  chapters  on  energy  as  driving  force  (Crickmay  1974;  Bloom 
1978;  Embleton  and  Whalley  1979). 

Odum  (1978a: 14)  presents  a  variety  of  energy  interpretations  of  land- 
scapes. He  states  that  landscape  is  "the  memory  and  the  bank  by  which  the 
flows  of  energy  develop  and  maintain  the  capital  against  depreciation,  and 
some  for  developing  other  values  that  go  out  in  exchange  for  inflows."  An 
outline  of  the  method  of  calculating  embodied  energy  in  landform  and  other 
structures  is  given  and  is  used  in  this  study.  The  method  is  applied  to  a 
Florida  cypress  wetland  and,  converting  energy  units  to  dollars,  a  value 
of  $3837/acre -yr  was  found  for  this  natural  land  type. 

Energetics  have  also  been  applied  to  landscapes  where  man  dominates 
over  natural  features  as  reviewed  above.     Energy  principles,  developed  by 
H.   T.   Odum  (1971,   1983),  have  been  used   for  regional  planning  (Odum  et  al . 
1975;  Odum  and  Peterson  1972;  Odum  et  al .   1979).     Considered  are  topics  of 


carrying  capacity,  investment  attraction,  predictions  of  future  scenarios, 
and  other  topics.     Applications  have  been  made  to  south  Florida   (Odum  and 
Brown  1979)  and  to  many  other  counties,  regions,  states,  and  countries. 
Costanza  (1975)  drew  maps  of  energy  inputs  (sun,  wind,   two  kinds  of  rain 
energy,  tides,  waves,  fossil  fuels,  and  electricity)  and  of  resultant 
structure  (biomass,  water  storage,  nutrients,   and  capital  structure  main- 
tained by  humans)   in  south  Florida  and  found  spatial  correlation  with 
economic  development.     The  structure  in  landform  was  not  evaluated, 
however,  which  might  have  improved  the  correlation  between  input  energy 
and  structure.     Using  the  energy  maps  as  a  base,   Costanza  (1979)  developed 
a  spatial  energy  exchange  model  that  generated  patterns  of  development 
that  approximated  real  landscape. 

In  other  energy  studies  of  human  landscape,  Woldenberg  (1968)  used 
the  maximum  entropy  concept  to  investigate  the  nature  of  central  place 
hierarchies.     He  attempts  to  explain  this  model  of  spatial  organization 
with  the  steady  state  principle  similar  to  the  work  of  Leopold  and  Lang- 
bein  (1962).     The  hexagonal  lattice  of  central  places  is  viewed  as  a 
least-work  lattice,  which  represents  the  state  of  maximum  possible 
entropy.     A  larger  work  (Woldenberg  1969)  relates  patterns  of  spatial 
structure  to  hierarchy  in  a  variety  of  system  including  cities,  rivers, 
trees,  etc. 

Odum  and  Kroodsma  (1976)  discuss  a  power  park  concept  as  a  future 
landscape.     Emphasis  is  placed  on  matching  high  energy— intens ive  land  uses 
with  larger  areas  of  low  energy— intensive  land  uses.     Energy  intensity  of 
land  use  along  with  other  design  concepts  have  been  discussed  by  the 
Southwest  Florida  Planning  Team  (1972). 
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Grillo  (1960)  attempts  to  generalize  geometric  form  into  two  classes 
a  cracking  pattern  as  an  expression  of  potential  energy  and  an  erosion 
pattern  as  an  expression  of  kinetic  energy.     Examples  are  taken  from  land- 
scape, art,  and  urban  design. 

Reclamation 

In  the  past,  around  the  turn  of  the  century,  reclamation  meant  bring- 
ing a  "wasteland,"  such  as  flood-prone  lands  and  deserts,  into  economic 
production.     Existence  of  naturally  produced  value  was  often  overlooked. 
At  that  time  there  was  a  National  Bureau  of  Reclamation  in  U.S.  govern- 
ment.    The  meaning  of  reclamation  is  different  now.     Reclamation  involves 
the  repair  of  a  disturbed  site,  returning  the  value  lost  in  disturbance. 
Reclamation  is  sometimes  defined  in  the  general  sense  of  replacing  any 
system  that  has  value  and  is  contrasted  with  restoration  that  is  taken  to 
mean  replacement  with  a  system  that  is  equivalent  to  the  original  system 
that  existed  before  disturbance.     No  such  distinction  is  made  in  this 
study  and  reclamation  and  restoration  were  used  interchangeably. 

Surface  mining  and  erosion  of  agricultural  land  are  two  of  the  more 
common  terrestrial  disturbances  requiring  reclamation.     Terrestrial  recla 
mation  often  consists  of  the  following:   land  grading  and  manipulation, 
soil  preparation,   fertilization,   irrigation,   erosion  control,   and  seeding 
and  planting.     The  task  is  essentially  one  of  establishing  an  ecosystem 
that  will  provide  a  cover  over  disturbed  areas   in  the  least  amount  of 
time.     Although  natural  succession  will  usually  restore  an  area,  reclama- 
tion is  invoked  to  speed  this   inherently  slow  process.     Costs  can  be  as 
high  as  shown  in  Table  1.     The  least  cost  proposed  is  $300/acre  for  a 
Florida  pasture  and  the  greatest  is  about  $300,000/acre   for  a  mangrove 
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forest  in  Puerto  Rico.    Where  succession  is  allowed  to  restore  the  land 
there  are  no  direct  costs.     Productivity  lost  during  disruption  is  often 
omitted  from  reclamation  accounts.     Often  the  costs  of  reclamation  appear 
to  be  exorbitant  and  their  economic  value  is  questioned   (Singer  1977). 
Partly  because  of  lack  of  knowledge  about  succession  on  disturbed  sites, 
use  of  succession  as  a  reclamation  alternative  is  uncertain,  and  costs 
involved  are  unknown.     In  this  dissertation,  objective  techniques  of 
energy  analysis  were  used  to  consider  reclamation  alternatives   in  phos- 
phate mining  in  Florida. 


Background  and  Concepts 

In  sections  that  follow,  concepts,  hypotheses,  and  definitions  are 
given  for  relating  energy  to  landscape  forms  and  storages.     The  energy 
measures  defined  here  were  evaluated  in  the  research  on  theories  of  land- 
form  development. 

Energy  Storage  Value  of  Landforms 

In  these  studies  storage  of  landforms  was  assigned  three  energy- 
related  values:  actual  potential  energy  content  given  in  Calorie  equiva- 
lents, embodied  energy  given  in  solar  equivalent  Calories,   and  energy 
effects  given  as  the  ratio  of  effect  to  feedback.     These  values  are  com- 
pared in  Figure  2.     The  actual  energy  stored  is  that   found  by  converting 
it  entirely  to  heat.     This  is  not  a  measure  of  work  or  effect  because 
energies  are  usually  used  in  systems  in  amplifying  interactions  with  othe 
energy  flows  rather  than  being  degraded  into  heat  without  interaction. 
The  same  actual  calories  of  different  types  of  energy,   such  as  solar, 
fossil  fuel,  and  electricity,  will  not  amplify  the  same  amount  of 
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Energy  Amplifier  Effect  is 
the  Ratio  of  Amplified  Production 
Flow  P  to  Actual  Calories  of  Feed- 
back Flow  F 


Actual  Potential 
Energy  Content 
(Actual  Stored 
Calories) 


Embodied  Energy  of  Storage 
is  Sum  of  Energy  Used  of 
Quality  S  to  Develop  Storage 
Over  Time 


Figure  2.     Model  showing  different  energy  values  of  storage. 
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productive  work.     This  is  another  way  of  saying  that  different  energy 
types  have  different  qualities. 

Embodied  energy  is  one  measure  of  energy  quality.     This  has  been 
defined  as  "the  energy  used  in  developing  energy  of  higher  quality"  (Odum 
and  Odum  1980).     It  is  the  amount  of  one  type  of  energy  used  to  make 
another  type.     Embodied  energy  storage  is  the  total  energy  equivalents  of 
one  type  used  over  the  time  needed  for  development  of  the  storage.  This 
meaning  is  similar  to  the  economic  notion  of  capital  cost. 

Energy  multiplier  effect  is  another  measure  of  the  quality  of  energy 
as  its  effect  or  "release"  value.     Ideally  this   is  the  amount  of  energy  it 
can  release  or  amplify  in  a  process  involving  interaction  with  other  ener- 
gy flows.     This  is  difficult  to  estimate  since  many  energies  feed  back  and 
forward  in  complex  webs  of  converging  and  diverging  flow,   and  individual 
effects  can  not  be  separated  easily.     It  might  be  expected  that  embodied 
energy  is  correlated  with  energy  effect  ("energy  cost"  =  "energy  effect"), 
but  such  relationships  have  not  been  tested.     In  this  dissertation, 
attempts  are  made  to  quantify  these  energy  measures  of  value  for  storage 
in  landforms  under  a  diversity  of  different  climatic  regimes  and  vegeta- 
tional  influences. 

Though  not  often  considered  as  such,  the  form  of  land  is  a  long-term 
storage.     Odum  (1975:147)   states,  "geomorphological   form  and  information 
as  measured  by  capital  energy  investment  becomes  as  much  a  state  variable 
as  biomass   in  biology  or  water  in  hydrology."     The  actual  storage  value  of 
landform  is  its  heat  content.     Geomorphic   forms  are  either  depressions, 
like  lakes  and  wetlands,  or  elevations,   like  mounds  and  hills.  Their 
actual  potential  energy  content  can  be  estimated  as  potential  energy 
stored  against  gravity,  either  to  fill   in  for  depression  forms  or  to  erode 
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away  for  elevated  forms.     Other  potential  energy  may  be  in  chemical  sub- 
stances of  the  storages. 

Embodied  energy  of  a  land  form  is  the  energy  used,   in  equivalents  of 
one  type,  multiplied  by  the  development  time.     This  approach  to  quantify- 
ing embodied  energy  assumes  that  a  production  process  is  competitive  so 
that   for  any  given  landform  development  time  has  been  minimized  through 
the  interaction  of  the  ecosystem  with  the  energy  sources.     In  this  sense 
natural  selection  may  operate  on  the  non-living  landform  directly  through 
the  living  ecosystem. 

Several  problems  are  involved  with  estimating  the  energy  in  landform, 
or  any  other  long-term  storage,   such  as  biomass,  which  is  made  up  of  cumu- 
lative additions  over  time.     Eight  algorithms  for  predicting  total  energy 
used  by  a  system  have  been  proposed  based  on  the  inflows  of  energies  (Odum 
1983).     No  one  method  is  apparently  most  appropriate.     In  this  study, 
energy  used  in  developing  landform  was  calculated  as  the  solar  energy 
embodied  in  wind,  rain,  river,  waves,  tides,  and  other  energies.  This 
approach  avoids  double  counting  of  sunlight   from  which  most  other  natural 
energies  are  derived  from  biospheric  cycles. 

Criteria  for  choosing  the  point   in  time  when  development  ends,  in 
order  to  calculate  embodied  energy,  are  not  clear.     Some  parameters  may 
continue  to  develop  after  others  have  stopped;   for  example,    in  a  forest 
the  diversity  of  species  may  continue  growing  after  biomass  has  reached 
steady  state.     Moreover,  does  energy  keep  becoming  embodied  even  after  a 
steady  state  is  reached  with  long-range  slow  changes  in  biological  evolu- 
tion?    For  the  purposes  of  calculations  here,   the  age  of  a  landform  was 
used  as  development  time  and,   for  ecosystems,   the  length  of  time  from 
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initiation  of  succession  until  climax  is  reached.     Climax  was  used  as  a 
period  when  main  storages  were  little  changing. 

Energy  effect  of  landform  is  the  amount  of  control  exerted  over  land- 
scape processes.     One  approach  to  quantifying  energy  effect   is  to  measure 
power  use  when  storages  are  being  utilized  without  replacement.     For  ele- 
vated land  forms  this   is   the  erosion  rate  and  for  depression  land  forms  it 
is  the  deposition  rate.     If  not  continually  renewed,  both  of  these  rates 
represent  loss  of  volume  of  storage.     The  landform  is  actually  degraded  by 
these  processes,  which  represents  work  that  yields  energy  effect. 

The  development  of  natural  vegetation  patterns  serves  to  control  ero- 
sion and  deposition  in  undisturbed  conditions.     Human  occupation  of  land- 
scapes, however,  generally  converts  natural  vegetation  to  other  uses, 
arrests  succession  or  substitutes  less  mature  vegetation  such  as  agricul- 
ture and  lawns.     The  effect  is  often  accelerated  erosion  in  one  area 
accompanied  by  deposition  elsewhere.     Man's  overall  impact  in  consumptive 
use  may  be  the  leveling  of  landscape  to  more  even  surfaces.     This  may  not 
be  the  best  use  of  landform  energy  unless  the  cycle  of  growth  and  develop- 
ment maximizes  power  in  the  long  run.     Odum  (1972)   cites  a  similar  example 
showing  growing  human  populations  drawing  on  large  ancient  geochemical 
storages  that  were  adapted   to  biosphere  cycles  as  a  step  to  start  a  new 
system  of  industrial  man. 

Hierarchy  in  Landscape 

In  relating  energy  to  area,  a  hierarchy  of  storages  develops  to  maxi- 
mize useful  work.     Landform  hierarchies  are  maintained  by  convergences 
over  large  areas,  not  just  by  energy  over  individual  structures.  Analo- 
gous with  food  chains  of  species  populations  within  ecosystems,  landscapes 
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may  be  considered  as  "food  chains"  of  land  forms  in  spatial  dimension. 
Dilute,  undifferentiated  expanses  capture  the  main  energy  flow  and  focus 
or  concentrate  it  into  structured  landscape  systems.  In  response  to  the 
concentrating  flows,  the  landform  may  redistribute  high  quality  services 
back  to  the  surrounding  system.  General  aspects  of  energy  distributions 
in  hierarchy  are  described  and  explained  by  Odum  (1983)  and  Brown  (1980). 
In  this  dissertation,  spatial  hierarchies  of  landform  are  considered. 

Energy  convergence  in  landscape  may  be  roughly  proportional   to  the 
ratio  of  areas  of  different  type.     Each  high  quality  landform  has  a  col- 
lection area  contributing  energy  to  it  through  water  flows  and  other  con- 
vergent pathways.     More  area  allows  more  development  of  the  hierarchy  to 
be  supported.     Recognizing  the  hierarchical  ratio  of  one  land  use  to 
another  helps  estimate  embodied  energies  where  pathways  converge.  For 
example,  Horton  (1945)  proposed  a  hierarchical  relationship  between  the 
relative  position  of  a  stream,   in  a  drainage  basin  (stream  order),  and  the 
log  of  number  of  streams  called  the  law  of  stream  numbers.     This  shows  a 
basin  with  many  small  streams  of  low  order  converging  energy  to  progres- 
sively fewer  streams  of  higher  order. 

Spatial  Expression  of  Energy  Sources 

Landscape  patterns  of  hierarchy  depend  on  the  spatial  distribution  of 
energy  sources.     Odum  (1983)  gives  three  types  similar  to  those  described 
below.     Also  a  fourth  is   included  here.     Figure  3  is  a  representation  of 
four  basic  spatial  distributions  that  energy  inputs  can  take — sheets, 
points,    fronts,  and  lines.     All  of  the  spatial  energy  distributions  act 
synoptically  over  relatively  large  areas  as  a  moving  force. 
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ENERGY  SHEET 


ENERGY  POINT 


ENERGY  FRONT 


ENERGY  LINE 


Figure  3.  Diagraramat ic  representations  of  four  basic  types  of  spatial  energy 
distributions.  Arrow  points  show  where  potential  energy  is  deliv- 
ered. 
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Energy  sheets  are  those  delivered  over  a  broad  area  perpendicular  to 
the  surface.     Sunlight  and  rainfall  are  examples  of  an  energy  sheet  from 
above  the  surface.     A  potentiometr ic  surface  of  groundwater  and  a  uniform 
rock  substrate  emerging  from  an  uplift   process  are  examples  of  a  sheet 
incoming  from  below  the  surface. 

An  energy  point  also  is  delivered  perpendicular  to  the  surface  but  at 
one  point  such  as  an  artesian  spring  rather  than  over  a  broad  area  as  in 
the  case  of  the  sheet. 

An  energy  front  is  defined  as  a  source  that  delivers  energy  along  an 
edge  that  moves  parallel  to  the  surface  and  sweeps  across  the  landscape. 
Examples  are  weather  fronts,  sheet   flows,   runoff,  and  tides.     These  ener- 
gies act  as  the  line  of  the  front  passes. 

Energy  lines  deliver  energy  along  a  line  parallel  to  the  land  sur- 
face.    Energy  is  supplied  by  force  along  the  line  and  at  the  end  point  of 
delivery.     Like  energy  points  the  line  force  is  expressed  over  a  small 
area.     Rivers  and  electric  power  lines  are  examples  of  energy  lines. 

Table  2  categorizes  many  natural  energies  into  the  classification 
scheme  of  spatial  expression  of  energy  sources. 

Defining  an  Energy  Source 

Each  energy  source  also  has  a  time  program  including  frequency, 
intensity,  duration,  and  direction.     To  fully  describe  an  energy  source 
then  requires  specification  of 

a.  actual  heat  equivalents  and/or  physical  units; 

b.  energy  type  and   its  embodied  energy  flow  expressed  in 
units  of  that  type; 


c.  spatial  distribution  of  energy  delivery;  and 


d.   program  of  delivery  with  time. 


Plan  of  Study 

Energy  basis  of  landforms  was  evaluated  including  those  that  are  part 
of  natural  landscapes  and  those  part  of  landscapes  of  human  economic 
development.     Quantitative  evaluations  of  embodied  energy  were  made  as  a 
measure  of  landscape  control  and  value.     An  energy  signature  of  interact- 
ing energy  inflows  was  calculated   for  ten  landforms.     For  each  a  minimodel 
was  constructed  of  the  landscape  structure  and  its  interaction  with  the 
input  energies.     Actual  energy  in  land  form  was  calculated  with  various 
potential  energy  equations.     Embodied  energy  in  storages  was  found  by  mul- 
tiplying the  energy  signature  used  over  the  support  area  by  the  energy 
transformation  ratios  and  the  development  time.     Embodied  energy  of  vege- 
tational  components  was  also  calculated  and  compared  with  land  form  values. 
Embodied  energy  of  landscape  was  related  to  location  of  human  developments 
and  to  aesthetics. 

The  theory  that  spatial  expressions  of  energy  sources  control  pat- 
terns of  landscape  was   tested   for   four  cases.     Landscape  pattern  was  used 
to  predict  dominant  energy,  which  was  then  calculated  independently.  A 
classification  of  the  shapes  of  landforms  is   presented,   and  the  adaptive 
value  of  shape  in  absorbing  energy  flows  is  discussed. 

The  energetics  of  developing  a  human-dominated  landscape  was  studied, 
using  the  phosphate  mining  district  of  central  Florida  as  an  example. 
Relative  values  of  human  and  natural   landscape  work  flows  were  compared. 
Energy  in  strip  raining  topography  was  evaluated  and  compared  with  natural 
landforms.     Ecosystem  development  on  spoil  mounds  was   investigated  by  mea- 
suring characteristics  of  successional  systems  on  mounds  of  different  age, 
ranging   from  1  to  50  years  since  mining.     Forest  structure  and  composition 
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on  the  oldest  spoil  mounds  were  compared  with  a  mature  undisturbed  forest 
to  assess  the  degree  to  which  succession  has  naturally  restored  develop- 
ment.   A  model  of  succession  was  constructed  for  southern  coastal  plains 
old  field  succession  and  validated  with  data  from  spoil  mound  succession. 
Simulations  tested  alternatives  for  accelerating  succession.     Energy  anal- 
ysis was  used  to  evaluate  mechanical  reclamation  versus  natural  succession 
as  alternatives.     A  model  of  a  typical  fertilizer  operation  was  con- 
structed to  place  reclamation  costs   in  perspective  of  whole   industry  bud- 
gets and  to  calculate  net  energy  of  phosphorus.     Recommendations  based  on 
fieldwork  and  model  simulation  were  made  for  effective  ecosystem  restora- 
tion on  disturbed  lands. 


Site  Description 

Phosphate  Mine  Study  Sites 

The  study  sites  for  fieldwork  on  successional  phosphate  mine  ecosys- 
tems were  located  in  southern  Polk  County  in  the  central  Florida  phosphate 
district  (Fig.  4a).     Eight  sites  were  studied  ranging  from  1  year  to  50 
years  in  age.     All  sites  but  one  were  located  at  the  Payne  Creek  mine  of 
the  Agrico  Chemical  Company  (Fig.   4b).     The  remaining  site  (50  year)  was 
about  16  kilometers  north  of  the  Payne  Creek  mine  in  the  Christina  resi- 
dential development  of  the  International  Minerals  and  Chemical  Development 
Corporation.     Study  site  ages  were  determined  by  a  review  of  mining  rec- 
ords on  file  at   the  Agrico  Chemical  Company  office  and  by  discussion  with 
mining  officials.     All  sites  were  abandoned  mine  areas  dominated  by  rows 
of  spoil  mounds  or  plateaus.     The  mounds  were   formed  of  overburden  mater- 
ials (sands)   stripped  and  deposited  by  the  operation  of  dredges  or  drag- 
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lines.    Mound  heights  ranged  from  3  to  12  meters.     Data  were  reported  by 
age  of  the  mound.     Where  two  sites  of  the  same  age  were  studied,  their 
data  were  averaged. 

The  50-year-old  sites  differ  significantly  in  geomorphology  from  the 
younger  sites.     At  these  older  sites  spoil  material  was  deposited  in 
blocks  as  large  plateaus  rather  than  in  the  more  recent  configuration  of 
spoil  mound  rows  with  intervening  troughs.     So  as   to  minimize  the  effect 
of  the  differences  in  original  spoil  shape,  only  the  slopes  of  all  areas 
were  studied. 

Sites  in  San  Felasco  Hammock 

In  order  to  assess  the  degree  to  which  natural  succession  restores 
ecosystem  structure,  the  mesic  hammock  forest  on  the  San  Felasco  Preserve 
was  sampled  and  compared  to  the  oldest  ecosystems  on  spoil  mounds.  San 
Felasco  is  a  mature  or  late  successional   forest  and  was  used  as  a  "climax 
analog"  to  compare  with  phosphate  succession.     The   forest  is   located  about 
13  kilometers  northwest  of  Gainesville  in  north-central  Florida  off  State 
Road  232.     It  is   in  a  protected  area  of  state  preserve  associated  with  a 
large  sink  hole,  the  Devil's  Millhopper.     Many  ecological  studies  have 
been  conducted  on  this   forest  including  aspects  of  vegetation  (Monk  1960), 
light  transmittance  (Ewel  1969),  mineral  cycling  (Ewel  et  al .   1975),  and 
metabolism  (Lugo  et  al.   1978).     The  sampling  site  used   in  this  study  was  a 
shallow  slope  with  moist  soil  and  many  evergreen  species.     Although  San 
Felasco  Hammock  is  approximately  200  kilometers  north  of  the  central  Flor- 
ida phosphate  region,   the  vegetation  is  to  some  extent  similar  and  compar- 
isons were  thought  to  be  useful. 


Landform  Sites  for  Embodied  Energy  Calculations 

Table  3  lists  the  landform  types  chosen,   their  location,   and  special 
adaptations.     Examples  were  picked  to  cover  a  variety  of  environments  with 
different  dominant   forcing  functions.     Comparisons  were  then  made  between 
arid-humid  types,  boreal-tropical  types,  and  across  multitudes  of  energy 
inputs.     Literature  was  surveyed   for  each  landform  type  and  a  particular 
location  was  chosen  that  had  sufficient  background  data  to  use  in  energy 
calculations . 
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CHAPTER  2 
METHODS 


Energy  Analysis 

Energy  Circuit  Language 

The  energy  circuit  language  (Odum  1967,   1971,   1976,   1983;  Odum  and 
Odum  1981)  was  used  to  model  land  systems.     This  is  a  symbolic  language 
for  representing  whole  systems  in  diagrams.     It  has  similarities  with 
electrical  circuit  diagrams  but  has  a  richer  content.     The  language  is  a 
kind  of  visual  mathematics  since  each  symbol  configuration  has  a 
mathematical  expression.     Diagrams  using  the  energy  symbols  have  the  same 
information  as  matrices  and  sets  of  equations  but  go  beyond  by  containing 
the  constraint  of  energy  laws  of  conservation  and  entropy  increase. 
Models  of  this  type  can  be  used  for  both  kinetics  and  energetics.     Use  of 
the  energy  circuit  language  forces  one  to  include  important  features, 
lessening  the  chance  that  something  significant  will  be  left  out  of  the 
model.     Artistic  use  of  hierarchy  is   incorporated   into  the  models  by 
arranging  the  position  of  symbols.     Units  having  low  energy  concentration 
(low  quality)  are  drawn  at  the  base  of  circuits  while  high  energy  units 
(high  quality)   are  drawn  at  the  end  of  circuits. 

The  general  procedure  in  diagramming  a  system  was  to  start  by  includ- 
ing the  main  energy  sources  then  adding  storage  and  pathways  along  which 
energy  flows   forward  and  back,   finally  leaving  the  system  as  dispersed 
heat  or  export  to  another  system.     Energy  was  used  as  a  common  denominator 
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since  all  real  flows  and  storages  have  energy  value  and  because  energy  is 
the  basis  of  causal  action. 

Expressing  Energy  in  Equivalent  Units 

After  diagramming,  energy  flows  and  storages  were  evaluated   for  their 
actual  energy  (heat)  content.     Equations  for  calculating  actual  energies 
are  given  below.     Variables  used  in  the  equations  are  listed  in  appendix 
tables  for  each  of  the  land  forms  analyzed.     Much  of  these  data  were  taken 
from  a  handbook  of  energy  values  for  the  United  States  (Odum  et  al .  1978). 

1 .  Solar  Insolation 

Values  reported  are  surface  measurements. 

2 .  Wind  Kinetic  Energy 

E  =  0.5pV2c(l/d)(2.38  x  10"11  Cal/erg) (3 . 15  x  107  s/yr) 

(104  Cm2/m2); 
p  =  density  of  air  =  1.2  x  10~3  g/cm3; 
V  =  wind  velocity  in  cm/ s ; 

C  =  eddy  diffusion  coefficient  =  104  cm2/s;  and 
d  =  height  of  boundary  layer  =  104  cm. 

3.  Precipitation  and  Runoff  Chemical  Potential  Energy 
E  =  nRT  ln(C2/C1)M; 

n  =  1  mole/18  g  H20 ; 

R  =  gas  constant  =  (1.99  cal/°mole)(l  Cal/1000  cal) 

=  1.99  x  10"3  Cal/°mole; 
T  =  absolute  temperature  in  °K; 
C2  =  water  content  of  precipitation  or  runoff 

=  1,000,000  ppm  -  TDS  of  precipitation  or  runoff; 
Cj  =  water  content  of  receiving  water 
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=  1,000,000  ppm  -  TDS  of  receiving  water; 
TDS  =  total  dissolved  solids  in  ppm;  and 
M  =  mass  of  precipitation  or  runoff  per  year  in  g/ra2*yr 
=  precipitation  in  cm/yr  (1  g/cm^)(10^  cm2/m2). 

4.  Precipitation  Impact  Kinetic  Energy 
E  =  0.5Mv2(2.38  x  10"11  Cal/erg); 

M  =  mass  of  precipitation  per  year  in  g/m2-yr 

=  precipitation  in  cm/yr  (1  g/cm-^)(10^  cm2/m2);  and 
v  =  velocity  of  raindrops  =  762  cm/s   for  an  average  drop  diameter 
(Odum  et  al.  1978). 

5.  Runoff  Elevated  Potential  Energy 
E  =MGHR(2.38  x  10"11  Cal/erg); 

M  =  mass  of  runoff  per  year  in  g/m2,yr 

=  runoff  in  cm/yr  (1  g/cm3)(104  cm2/m2); 
G  =  gravitational  acceleration  =  980  cm/s2;  and 
HR  =  elevation  drop  or  head   in  cm. 

6.  Tidal  Potential  Energy 

E  =  0.5pGH^n(2.38  x  10"11  Cal/erg)(104  cm2/m2); 

p  =  density  of  seawater  =  1.025  g/cm3 ; 

G  -  gravitational  acceleration  =  980  cm/s2; 

Hj  =  tidal  head  in  cm;  and 

n  =  number  of  tides  per  year. 

7.  Organic  Inflow  Potential  Energy 
E  =M0c; 

MQ  =  mass  of  organic  matter  per  year  in  g/m2-yr; 

c  =  heat  content  of  organic  matter  =  4.5  Cal/g  (E.   P.   Odum  1971). 
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8.  Wave  Potential  Energy 

E  =  l/8pG3/2H^/2w(2.38  x  10-11  Cal/erg)   (3.15  x  107  s/yr) ; 

p  =  density  of  seawater  =  1.025  g/cm3; 

G  =  gravitational  acceleration  =  980  cm/s2; 

Hw  ■  wave  height  in  cm;  and 

W  =  width  of  area  affected  in  m. 

9.  Flood  Head  Potential  Energy 

E  =  pDGHFF(l/A)(2.38  x  10"11  Cal/erg); 

p  =  density  of  water  =  1  g/cm3; 

D  =  flood  discharge  in  cm3/s; 

G  =  gravitational  acceleration  =  980  cm/s2; 

Hp  =  flood  head  or  height  in  cm; 

F  =  duration  of  flood  per  year  in  s/yr;  and 

A  =  area  of  flooded  zone  in  m2. 

10.  Uplift  Potential  Energy 

E  =  MGHU  (2.38  x  10-11  Cal/erg); 

M  =  mass  of  uplifted  material  per  year  in  g/m2-yr 

=  height  of  uplift  in  cm/yr  (  ,  density  of  material) (10^  cm2/m2) 
G  =  gravitational  acceleration  =  980  cm/s2;  and 
Hu  =  height  of  uplift  per  year  in  cm. 

1 1.  Storm  Kinetic  Energy 
E  =  Sn; 

S  =  kinetic  energy  per  storm  in  Cal/m2;  and 
n  =  number  of  storms  per  year. 

12.  Freeze-Thaw  Cycle  Energy 
E  ~  Hp_*j*n ; 

Hp_T  =  heat  energy  per  freeze-thaw  cycle 


=  [(AT)(HC)(VF)  +  (Hm)(VF)](fiT/T)/At; 
AT  =  temperature  change  over  freeze-thaw  cycle  in  °C ; 
T  =  absolute  temperature  in  °K; 

Hc  -  heat  capacity  of  soil  substrate  in  Cal/cm3-°C; 

1^  =  heat  of  melting  =  8.0  x  10"2  Cal/cm3; 

Vp  *  volume  affected  by  freeze-thaw  cycle  in  cm3/m2; 

At  =  length  of  cycle  in  days;  and 

n  ■  number  of  daily  cycles  per  year. 

13.  Sediment  (sand  or  mixture)  Potential  Energy 
E  =  MGH; 

M  =  mass  of  accumulated  sediments  per  year  in  g/m2*yr 

=  PL (104  cm2/m2); 
p  =  density  of  sediments  in  g/cm3; 
L  =  accumulation  rate  in  cm/yr; 
G  =  gravitational  acceleration  =  980  cm/s2;  and 
H  =  0.5L. 

14.  Actual  Potential  Energy  in  Landforms 

E  =MGH(2.38  x  10_11  Cal/erg)(104  cm2/m2); 
M  =  mass  of  land  form  =  pV; 

p  =  density  of  landform  materials  in  g/cm3; 
V  =  volume  of  a  specified  shape  in  cm3; 
G  =  gravitational  acceleration  =  980  cm/s2;  and 
H  =  average  height  of  landform  in  cm. 

To  convert  energies  to  a  common  basis  of  equivalent  work  units 
embodied  energy,   actual  energy  values  were  multiplied  by  energy 


transformation  ratios  (quality  factors)  that  are  ratios  of  Calories  of 
solar  energy  per  Calorie  of  another  type. 

Energy  transformation  ratios  are  energy  transformation  efficiencies 
defined  as  the  amount  of  energy  of  the  type  required  to  develop  unit 
energy  of  another  type.     Energy  transformation  ratios  used  in  this  study 
are  given  in  Table  4,   primarily   from  Odum  and  Odum  (1980:7),  who  noted, 
"at  this  stage  in  energy  analysis  calculations  are  made  objectively  and 
documented  with  the  realization  that   tables  of  conversion  efficiencies 
will  be  improved  with  time."     Other  energy  transformation  ratios  were 
calculated  from  case  studies.     Solar  energy  was  used  as  the  base  quality 
unit  in  all  calculations. 

An  overall  ratio  between  energy  and  money  was  used  to  estimate  energy 
embodied   in  high-quality  goods  and  services.     It  was   found  by  dividing  the 
gross  national  product  into  the  total  energy  use  by  the  general  economy, 
including  purchased  fuels  and  natural  renewable  energy,  both  in  energy 
equivalents  of  the  same  quality.     This  ratio  is  the  amount  of  energy  used 
per  unit  dollar  circulating  in  the  economy,   which  takes   into  account  the 
free  natural  energies  contributing  to  the  economy.     For  the  United  States 
in  1980,   this  ratio  was  11,000  coal  equivalent  Calories   per  dollar  (Odum 
and  Odum  1980). 

Parameters  Used  for  Energy  Analysis 

From  an  elevated  model,  analysis  was  made  of  the  magnitudes,  distri- 
butions, and  proportions  of  storages  and   flows  of  energies.     All  energies 
were  translated  into  embodied  energy  equivalents  of  one  type.     Values  and 
indices  were  calculated  from  the  evaluated  model  to  aid  in  assessment  of 
sources,   storages,   and  alternative  actions  as  shown  in  Figure  5. 
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Table  4.     Energy  transformation  ratios  used  to  convert  different  energy 
sources  into  equivalent  units  of  solar  energy  (from  Odum  and 
Odum  1980;  Odum  et  al .  1978). 


Transformation  Ratios, 
Energy  Source  Direct  Solar  Cal/Cal 


Sun 

1 

12.6 

Wind 

Thunderstorms3 

1100 

Coal 

2000 

Organic  inflows 

2400 

Waves 

2700 

Elevated  potential  of  streams 

3000 

Petroleum  and  oil 

3100 

Tide,  physical  potential 

3200 

Chemical  potential  of  rain 
over  land 

3600 

Electricity 

8000 

Kinetic  impact  of  rain 

2.4 

x  105 

Flood  currents 

3.6 

x  105 

Geologic  uplift 

4.4 

x  1012 

'Calculated  in  this  study. 
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S=  SOLAR  ENERGY  USED 
A|  AzAz  =  AUXILIARY  ENERGIES 

E=  OUTSIDE  CONTROL  ENERGY 
0=  OUTPUT  FOR  SALES  OR 
EXCHANGE 

Figure  5.     Definitions  used  for  energy  analysis.     Energy  signature  =  S,  Ax , 
A2   •    •    •  An'  a11   in  units  of  the  same  quality.     Net  energy  of  0 
=  0  -  E,  both  in  units  of  the  same  quality.     Yield  ration  of  0  = 
0/E.     Total  embodied  energy  flow,  Ec  =  S  +  [   (AM  -  S)   if  posi- 
tive]  +  An.     Where  A^j  is  energy  of  indirect  solar  origin  and 
An  is  of  non-solar  origin  (such  as  tide). 
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The  energy  signature  is  the  set  of  all  energy  sources  including 
solar,  which  is  the  principal  natural  external  source,  auxiliary  energies 
such  as  wind,  rain,  and  uplift — which  are  by-products  of  solar — and  pur- 
chased energy  such  as  goods,  fuels,   services,  and  capital  investment  from 
human  economies. 

The  solar  energy  portion  of  the  energy  signature  was  calculated  as 
the  direct  solar  input  plus  the  indirect  (auxiliary)   solar  energy  of  the 
oceans  and  atmosphere  driving  weather  and  geological  energies.     The  direct 
solar  was  first  subtracted  from  each  indirect  energy  to  eliminate  double 
counting  and  any  positive  remainder  was  added  back  to  solar.     This  is  the 
energy  available  to  be  used  in  landscape  work.     The  sum  of  direct  and  cor- 
rected indirect  energies  was  used  to  avoid  counting  auxiliary  energy  twice 
since  they  are  by-products  of  direct  sunlight. 

Embodied  energy  in  the  long-term  storage  of  landform  was  calculated 
by  the  equation  given  below 

Landform  embodied  energy  =  ETCDT  (2.1) 

where  ET  =  total  embodied  energy  input  in  solar  Cal/m2-yr;  C  =  concentra- 
tion area  in  m2;  and  DT  =  landform  development  time  in  years.  Values 
for  the  total  embodied  energy  input,  ET,  are  given  in  tables  summarizing 
the  energy  signatures   for  each  landform.     Data  on  concentration  area,  C, 
and  development  time,  Dj,  are  given  in  appendix  tables   for  each  landform. 

Net  energy  of  an  output  is  the  difference  between  the  embodied  energy 
of  the  output  and  the  exchange  energies  used  to  extract  the  output.  The 
yield  ratio  of  this  extraction  process  is  the  output  energy  divided  by  the 
exchange  energies. 
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Successional  Model 

A  model  of  succession  was  simulated  to  test  alternatives  for  acceler- 
ating ecological  development.    Modeling  was  done  using  the  mathematical 
expressions  that  accompany  the  energy  circuit  language.     Figure  6  shows 
the  model  along  with  the  expressions  for  storages  and  pathways.  Three 
stages  of  succession  are  shown  with  storages  corresponding  with  biomass: 
Q3  is  grasses,  Q2  is  pines,  and  Q:  is  hardwoods.     Q4  is  litter 
storage  that  receives  input  from  all  stages  of  succession. 

The  main  driving  energy  source  was  sunlight,  which  is  apportioned 
according  to  the  maximum  stature  and  efficiency  of  each  successional 
stage.     Hardwoods  capture  the  most  sunlight  since  they  grow  to  be  the 
tallest  and  are  the  most  efficient  with  only  the  sunlight  remaining  being 
available  to  pines  and  grasses.     Pines  get  the  next  use  of  sunlight  fol- 
lowed by  the  grasses,  which  get  the  remainder   from  the  pines.     This  was  a 
way  of  simulating  shading  effects  that  incorporate  kinetics  of  limiting 
availability  of  sunlight. 

Initially  storages  were  started  at  zero  and  each  received  germination 
inputs   from  outside  seed  sources.     Seeding   input   to  pines  and  hardwoods 
were  controlled  by  threshold  switches  of  the  litter  storage.  These 
switching  controls  on  germination  were  related  to  changing  microclimates 
of  the  forest   floor,  which  allow  germination  of  seeds.     An  additional 
threshold  of  light  level  was  also  incorporated  as  control  on  hardwood  ger- 
mination.    The  model  was  run  with  and  without  this  switch  for  comparison 
when  light  was  potentially  limiting  to  germination. 

The  model  was  calibrated  with  data  largely  from  Georgia  old  field 
studies.     Timing  of  stages  was  adjusted  through  germination  thresholds  to 
Florida  phosphate  succession,  which  is  much  faster  than  Georgia  old  field 
succession.     The  model  was  simulated  with  a  1 -quarter-year  (3  month)  time 
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interval.     Smaller  time  steps  were  used  in  runs  simulating  litter  data  for 
validation.     Simulations  were  run  on  an  Intecolor  8000  desk-top  computer 
using  BASIC. 

Ecological  Field  Measurements 

Ecosystem  Parameters 

Percent  vegetation  cover  was  estimated  visually  from  0.25-m2  plots. 
Cover  measurements  were  taken  every  3  meters  over  transects  ranging  in 
length  from  6  to  36  meters  depending  on  the  particular  conditions  of  the 
surface.     Sample  sizes  were  60  plots  on  1-year-old  mounds,  92  plots  on  5- 
year-old  mounds,  68  plots  on  10-year-old  mounds,   and  84  plots  on  20-year- 
old  mounds.     On  the  oldest  sites  cover  was  estimated  from  aerial  photo- 
graphs.    Optical  density  as  an  index  of  biomass  (Odum  et  al.   1963a)  was 
measured  with  a  portable  solar  radiometer  (Matrix  Inc.,  Mark  IV  solarim- 
eter),  which  reads  percent  transmittance  between  a  setting  made   in  a 
clearing  and  that  under  canopy.     Sample  sizes  were  73  on  1-year-old 
mounds,   73  on  5-year-old  mounds,   59  on  10-year-old  mounds,   80  on  20-year- 
old  mounds,  84  on  50-year-old  mounds,  and  60  at  San  Felasco.  Percent 
transmittance  was  converted  to  optical  density  by  dividing  by  the  cosine 
of  the  angle  between  the  position  of  the  sun  and  the  vertical  position, 
then  taking  the  log  and  changing  sign.     Because  of  the  trigonometric 
transformation  no  estimates  of  embodied  energy  were  calculated.  Measure- 
ments were  taken  under  vegetation  at  ground  level.     Tree  height  was  mea- 
sured for  canopy-level   trees  with  a  Haga  altimeter.     Sample  sizes  were  16 
trees  on  10-  and  20-year-old  mounds  and  20  trees  on  50-year-old  mounds  and 
at  San  Felasco.     Ten  litter  samples  were  taken  by  hand  from  0.25-m2 
plots  at  each  age  mound.     Samples  were  dried  at  60 °C   for  24  hours  and 
weighed.     Litter  depths  were  measured  whenever  litter  was  greater  than  0.3 
cm  deep.     At   lesser  depth,    litter  was   too  sparse   to  form  a  distinct  layer. 
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Sample  sizes  were  32  on  5-year-old  mounds,  20  on  10-yearold  mounds,  32  on 
20-year-old  mounds,   60  on  50-year-old  mounds,  and  64  at  San  Felasco. 

Soil  organic  matter  was  measured  by  loss  on  ignition.     Ten  soil  sam- 
ples were  taken  from  the  upper  10  cm  for  each  age  mound.     Subsamples  of 
these  were  dried  at  60°C  for  24  hours  and  weighed,  then  ashed  at  500°C  for 
1  hour  and  weighed  again.     Organic  matter  was  estimated  as  the  difference 
between  the  dry  weight  and  ashed  weight.     Organic  matter  values  obtained 
with  this  technique  may  be  overestimated  because  of  combustion  of  water 
bound  to  clays  that  is  lost   in  addition  to  combustion  of  organic  matter. 

Soil  profiles  were  also  described  from  the  youngest  (1  year  old)  and 
from  one  of  the  oldest  (50  year  old)   sites.     Descriptions  were  based  on 
two  one-meter-deep  and  two  half-meter-deep  soil  pits  at  each  site.  Ant 
nests  as  an  index  of  soil  animal  activity  were  counted  in  0.25-m2  plots 
used  for  vegetation  cover  determination.     No  estimates  of  standard  error 
were  calculated  since  so  few  plots  had  nests.     Data  are  reported  as  nests 
per  10  m2.     Sample  sizes  were  60  plots   (15  m2)  on  1-year-old  mounds, 
92  plots  (23  m2)  on  5-year-old  mounds,  and  87  plots  (21.75  m2)   on  10- 
year-old  mounds.     Because  of  the  density  of  vegetation  it  was  not  possible 
to  count  nests  at  sites  older  than  10  years. 

Trees  greater  than  10  centimeters  dbh  were  sampled   for  density,  and 
basal  area  was  calculated  at   the  50-year-old  sites  and  at  San  Felasco. 
Three  900-ra2  plots  at  the  site  in  the  Payne  Creek  mine  and  at  San  Fel- 
asco and  five  450-m2  plots  at  Christina  were  sampled.     The  smaller  plot 
size  at  Christina  was  due  to  the  shorter  length  of  the  slopes  (about  15 
meters  at  Christina  and  about  30  meters  or  more  at   the  Payne  Creek  mine 
and  at  San  Felasco).     Importance  values  were  used  to  assess  the  biological 
contribution  of  a  tree  species  to  the  forest  community.     Importance  was 
defined  as  the  sum  of  relative  density  plus  relative  dominance  (as  mea- 
sured by  basal  area  in  m2/ha)  divided  by  two  to  put  the   index  on  a  per- 
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centage  basis.     This  definition  omits   frequency  that  was  originally  in- 
cluded by  Curtis  and  Mcintosh  (1951).     Because  of  the  small  number  of 
plots  and  their  large  size,  frequency  values  were  high  and  would  inflate 
importance  percentages  (see  also  Daubenmire  1968).     A  measure  of  relative 
frequency  in  terms  of  percent  occurrence  in  the  total  number  of  plots, 
however,  was  included  separately. 

Indices  for  Comparison  of  Forests 

One  of  the  goals  of  the  field  studies  was  to  determine  the  extent  to 
which  succession  has  restored  the  natural  ecosystems  on  spoil  mounds.  One 
way  this  was  accomplished  was  by  comparing  indices  of  tree  community 
structure  of  the  oldest  phosphate  forests  with  a  "climax  analog,"  San  Fel- 
asco  Hammock,   near  Gainesville.     Five  indices  were  used — the  Holdridge 
Complexity  Index  (Holdridge  1967),   species  per  1000  individuals  (Odum, 
Cantlon,   and  Kornicker  1960;  Margalef  1958),   the  Shannon  Index  (Shannon 
and  Weaver  1963;  Margalef  1958),  and  the  Simpson  Index  (S  impson  1949) 
along  with  a  floristic  similarity  index  (Sorensen  1948). 

The  Holdridge  complexity  index  is  calculated  as  follows: 

HBDS 

CI=1000  (2-2) 

where  H  =  average  height  of  the  three  tallest   trees  (m);  B  =  basal  area 
(m2/ha);  D  =  density  (number  of  stems  of  10  centimeters  dbh  or  great- 
er/ha); and  S  -  number  of  species  (of  10  centimeters  dbh  or  greater). 
These  values  are  usually  extrapolated  from  a  0.1-hectare  sample  area,  but 
in  this  study  data  were  based  on  the  total  sample  areas  of  0.5  hectares  at 
the  phosphate  forest  and  0.27  hectares  at  San  Felasco. 

Cumulative  species  per  1000  individuals   is  an  expression  of  the 
approximate  logarithmic  distribution  of  commonness  and  rarity  in  natural 
communities.     The  index  is  usually   found  by  plotting  the  cumulative 
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appearance  of  new  species  as  a  function  of  the  logarithm  of  the  cumulative 
number  of  individuals  examined.     Here  species  per  1000  was   found  by  back 
extrapolating  from  a  straight  line  plot  passing  through  the  origin  and  the 
point  established  by  the  log  of  the  total   individuals  counted  versus  the 
number  of  species  found. 

The  Shannon  Index,  derived   for  the  theory  of  communication,  has  often 
been  used  as  a  measure  of  diversity  in  ecological  communities.  Its 
expression  is 

H  =  -2>i>log2(Pi>  (2.3) 
where  Pj_  ■  (nj;)/N  =  observed  frequency  for  each  species;  n^  =  number 
of  individuals  for  each  species;  and  N  =  total  of  all   individuals.  The 
units  of  this  index  calculated  with  log2  are  in  "bits  per  individual." 

Simpson's   Index  is  defined  as 

C  =  1  -£(ni/N)2  (2.4) 
with        and  N  as  above.     This  index  ranges  from  a  minimum  diversity  of 
zero  with  one  individual  of  one  species   toward  a  maximum  diversity  of  one 
when  n£  is  small  and  N  is  large. 

Similarity  between  floristic  species  composition  was  measured  with 
the  Sorensen  Index  given  below 

S  =  [2w/(a  +  b)]   x  100  (2.5) 
where  a  =  number  of  species  in  community  a;  b  =  number  of  species  in 
community  b;  and  w  =  number  of  common  species  between  communities  a  and  b. 
This  index  has  the  form  of  a  percentage  ranging   from  a  minimum  of  zero 
when  there  are  no  common  species  (w  =  0)   to  100  when  all  species  are  held 
in  common  (w  =  a  =  b) . 


CHAPTER  3 
RESULTS 


Embodied  Energy  of  Natural  Landforms 

For  each  type  of  land  form  developed   in  geological,  ecological,  and 
cultural  work,  a  model  is  given  showing  the  main  driving  energies,  system 
processes,  and  storages.     The  energy  signature  is  shown  with  a  table  of 
embodied  energy  flows.     Then  data  are  given  for  representative  storage 
parameters  and  development  times  involved.     The  role  of  these  data  tables 
was  to  document  the  range  of  variability  of  landform  dimensions  and 
process  rates  in  reference  to  the  site  chosen  for  case  study.  Finally 
results  of  calculations  estimating  magnitudes  of  embodied  energy  in 
resulting  landforms  are  given. 

Embodied  Energy  in  Floodplain 

A  floodplain  is  a  landform  organized  by  river  flooding.     Many  flood- 
plains  are  composed  of  sediments  deposited  by  river  floods.     These  have 
much  organic  sediment  from  floodplain  vegetation.     Vegetation  of  flood- 
plains  is  usually  better  developed  than  its  surroundings  due  to  rich  soils 
and  a  longer  hydroperiod.     A  model  of  a  floodplain  is  shown  in  Figure  7. 
The  flood  is  shown  as  a  switch  operated  by  a  comparison  of  water  volume 
and  channel  volume.     When  water  volume  is  greater  than  channel  volume, 
flooding  occurs  leading  to  deposition  of  sediments,  which  make  up  the 
floodplain  landform.     Table  5  provides  background  data  on  floodplains. 
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Average  deposition  is  about  1  foot  per  major  flood.     For  case  study  the 
Huron  River  in  southeastern  Michigan  was  chosen.     The  energy  signature  of 
the  Huron  River  floodplain  is  given  in  Table  6.     The  dominant  energies  are 
organic  materials  deposited  by  floods  and  the  physical  energy  in  the  flood 
itself.    Using  a  development  time  of  10,000  years  (time  since  glacier 
recession)  gives  the  embodied  energy  of  1.0  x  10^  solar  Cal/m^. 

Embodied  Energy  in  Coral  Reef 

Coral  reefs  are  found  throughout  shallow  tropical  water  that  is  low 
in  nutrients.    The  system  is  mostly  nonliving  calcium  carbonate  skeletons 
of  corals  and  algae  with  a  covering  of  living  biomass.     A  coral  reef  model 
is  shown  in  Figure  8.     The  main  storage  is  reef  form,  which  interacts  with 
waves  and  tides  to  form  turbulence.     Because  of  the  concentrated  biomass, 
turbulence  allows  the  exchange  of  oxygen,  disposal  of  wastes,  and  acquisi- 
tion of  food  and  nutrients   from  the  surrounding  waters.     Calcium  input  is 
shown  but  was  not  evaluated.     Background  data  on  reef  growth  and  current 
velocity  for  several  reefs  are  given  in  Table  7.     An  energy  signature  for 
the  Florida  Key  reefs  is  given  in  Table  8  and  was  used  as  a  case  study. 
For  a  7-meter  section,  as  in  the  Pleistocene  Florida  Key  reef,  embodied 
energy  was  found  to  be  5 .  1  x  10^  solar  Cal/m^. 

Embodied  Energy  in  Salt  Marsh 

Along  low  gradient  sedimentary  coasts,  marshes  develop  with  salt- 
adapted  grasses.     In  certain  areas  of  resistant  rock,  marshes  can  develop 
on  their  own  peat  with  little  inorganic   sediment.     A  salt  marsh  model  is 
shown  in  Figure  9.     Turbulence  in  marsh  waters  from  tide  causes  transport 
of  inorganic  sediments  and  organic  matter  out  of  marsh.  Vegetation 
decreases  turbulence  by  friction  and  leads  to  deposition  of  sediment  that 
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Table  7.     Coral  reef  characteristics. 


Location 

Reef  Growth, 
cra/yr 

Current 
Velocity , 
cm/  s 

Re  ference 

Range  for  Caribbean  Reefs 

0.9-1.5 



Adey  1978 

Eniwetok 

1.6a 

0.49 

Odum  and  Odum 
1955 

Pleistocene  Reef,  Florida 
Keys 

0.43-0.85 

Ho  f  fme  i  s  t  er  and 
Multer  1964 

Grand  Cayman  Islands 

0.35 

Roberts  et  al . 
1976 

aMore  recent  evidence  suggests  that  no  net  reef  growth  may  have  occurred 
over  the  last  few  thousand  years  (Smith  1973). 
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comprise  the  form  of  the  marsh.     Table  9  gives  background  data  on  salt 
marshes.     Accumulation  rates  range  from  0.2  cm/yr  to  20  cm/yr,  indicating 
a  range  of  energy  conditions.     Using  Sapelo  Island,  Georgia,  as  a  case 
study,  an  energy  signature  of  a  marsh  is  given  in  Table  10.     Tidal  ener- 
gies are  highest,  providing  flushing  actions  and  building  intricate  chan- 
nel networks.     Assuming  a  development  time  of  300  years,   the  embodied 
energy  of  the  salt  marsh  is  9.7  x  109  solar  Cal/m2. 

Embodied  Energy  in  Mountain  Glacier 

Glaciers  are  massive  storages  of  ice  found  in  polar  and  alpine  re- 
gions.    A  model  of  a  mountain  glacier  is  shown  in  Figure  10.     Main  stor- 
ages of  the  glacier  are  snow  and  ice.     The  transformation  of  snow  to  ice 
is  shown  as  a  function  of  pressure.     The  important  role  of  glacier  ice  in 
erosion  is  shown  as  a  drain  on  mountain  storage.     However,   some  of  this 
energy  cycles  back  to  amplify  uplift  and  mountain  building.  Background 
data  on  glaciers  are  given  in  Table  11.     Blue  Glacier  in  the  Cascade  Moun- 
tains of  Washington  State  was  chosen  as  a  case  study.     A  general  energy 
signature  for  the  Cascades  is  given  in  Table  12.     Uplift  is  the  largest 
energy  source.     Using  a  development  time  of  30,000  years,  the  embodied 
energy  of  the  Blue  Glacier  was  calculated  at  2.6  x  lO1^  solar 
Cal/m2. 

Embodied  Energy  in  Palsa 

Palsa  are  ice-cored  mounds  found  in  northern  bogs  with  discontinuous 
permafrost.     They  have  a  thick  cover  of  moss  and  peat,  which  insulates  and 
maintains  the  mound.     Ice  is  concentrated   from  the  surrounding  bog  waters 
and  its  expansion  causes  the  mound  to  be  elevated.     Figure  11  shows  a 
model  of  a  palsa.     Vegetation  and  peat  are  shown  as  a  drag  on  heat  inputs 
and  outputs  from  the  mound.     Flows  between  water  and   ice  in  the  mound  are 
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shown  as  controlled  by  heat  storage.     Ice  grows  with  heat  loss  to  the 
atmosphere.     Table  13  provides  background  data  on  palsa  height  and  age. 
Much  controversy  exists  concerning  the  formation  of  palsa  as  indicated  by 
the  wide  range  in  ages.     For  case  study,  data  for  Devon  Island  were  used. 
An  energy  signature  of  Devon  Island  is  given  in  Table  14.     As  has  been 
pointed  out   for  these  types  of  systems   (Sigafoos  1952),   the  dominant 
energy  is  in  the  freezing  and  thawing  of  soils  and  waters,  which  lead  to 
the  concentration  of  ice  in  the  mound   form.     Embodied  energy  in  the  Devon 
Island  palsa  is  9.2  x  109  solar  Cal/m2. 

Embodied  Energy  in  Sand  Dune 

In  deserts  and  along  coasts  with  ample  sand  supply  and  little  vegeta- 
tion, sand  accumulates  in  sculptured  mounds  or  dunes.     Some  dunes  are  fix- 
ed, while  others  move  by  cumulative  saltation  across  landscapes.     A  model 
of  a  sand  dune  is  given  in  Figure  12.     Sand,  carried  by  wind,  is  deposited 
to  form  the  dune  by  the  action  of  vegetation  slowing  wind  velocity.  Back- 
ground data  are  given  in  Table  15.     Accretion  rates  range  from  3  to  200 
cm/yr.     Using  White  Sands,  New  Mexico,  as  a  case  study,  an  energy  signa- 
ture was  calculated  (Table  16).     Although  dunes  are  often  thought  to  be 
built  by  wind  energy,   sand  is   found  to  be  the  largest  embodied  energy 
input.     Embodied  energy  of  the  White  Sands  National  Monument  dunes  was 
found  to  be  3.0  x  1010  solar  Cal/m2. 

Embodied  Energy  in  Arroyo 

In  the  southwestern  regions  of  the  United  States  and  to  some  extent 
in  other  arid  regions,  sharply  cut  drainage  channels  called  arroyos 
develop.     Their  origin  has  been  related  to  sparse  vegetation  and  desert 
rain  regimes;  though  there  exists  some  controversy  as  to  the  actual  cause. 
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This  region  has  high  grazing  intensities  that  alter  natural  vegetation 
that  may  have  controlled  erosion  in  the  past.     A  model  of  an  arroyo  is 
given  in  Figure  13.     The  main  storage  is  channel  volume,  which  drains  the 
watershed  of  the  small  amounts  of  rain  characteristic  of  these  areas. 
Biomass  is  small  and  appears  to  exert  little  control.     Background  data  are 
given  in  Table  17.     The  extremely  fast  erosion  rates  are  notable, 
especially  considering  volumes  of  the  channels.    An  energy  signature  was 
calculated  (Table  18)   for  the  area  of  Rio  Puerco,  New  Mexico,  as  case 
study.     The  highest  energy  is  from  storm  rains.     Embodied  energy  of  the 
arroyo  was  calculated  at  5.1  x  109  solar  Cal/m2. 

Embodied  Energy  in  Oyster  Reef 

Oyster  reefs  are  dense  accumulations  of  oyster  shell   formed  in  turbid 
coastal  river  mouths  and  inter tidal  zones.     Unlike  coral  reefs,  these  are 
consumer  reefs  of  low  diversity.     They  are  adapted  to  somewhat  narrow 
salinity  ranges  of  about  25%o  to  15%o.     An  oyster  reef  model  is  shown 
in  Figure  14.     The  storage  of  shell  interacts  with  tide  forming  currents 
and  turbulence.     Currents  are  used  to  bring  in  food  and  remove  wastes,  as 
with  the  coral  reef.     Oysters   filter  plankton  and  detritus   from  surround- 
ing water  and  return  nutrients  by  way  of  recycle.     Background  data  for 
several  reef  parameters  are  given  in  Table  19.     Crystal  River,  Florida, 
was  used  as  a  case  study  with  some  data  from  other  locations.     An  energy 
signature  for  Crystal  River  is  given  in  Table  20.     Highest  energy  is  in 
the  plankton  organic  matter  followed  by  chemical  potential   in  runoff.  The 
latter  is  significant   in  that  it  represents  the  energy  necessary  to  con- 
trol the  salinity  range  in  which  oysters  are  found.     Using  the  Mobile  Bay 
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Table  19.     Oyster  reef  characteristics. 


Locat ion 


Intertidal  Crystal  River 
Mobile  Bay 

Texas 

Intertidal  North  Inlet 
Intertidal  Sapelo  Island 


Biomass , 

g/tn2  Age,  yr 


269.  1 


190-234 


310-500 
1108 


Thick- 
ness,  ft  Reference 


920 


Lehman  1974 

May  1971 
Bernard  1974 

Norris  1953 

Dame  1976 

Bahr  1976 
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reef  development  time,  the  embodied  energy  was  found  to  be  1.1  x  10^ 
solar  Cal/m^. 

Embodied  Energy  in  Florida  Bay  Mud  Mounds 

In  the  shallow  waters  of  Florida  Bay  off  the  southern  tip  of  peninsu- 
lar Florida  is  found  a  unique  type  of  carbonate  mud  mound.     The  mud  is 
largely  produced  by  carbonate-depositing  algae,  and  the  mounds  are  bound 
with  the  roots  of  seagrasses  on  which  the  algae  grow.     A  model  of  the 
Florida  Bay  mounds  is  shown  in  Figure  15.     Turbulence   from  wind,   tide,  and 
currents  tends  to  suspend  the  calcareous  mud  in  the  water  column.  How- 
ever, the  seagrass  acts  to  decrease  turbulence  by  friction,  which  causes 
turbidity  to  settle.     Calcium  carbonate  particles  and  mound  form  are  built 
directly  by  algal  biomass.     Data  were  taken  from  Scholl   (1966)   and  Enos 
and  Perkins  (1979).     An  energy  signature  for  Florida  Bay  is  given  in  Table 
21.     The  highest  energy  is   from  chemical  runoff  from  drainage  of  the  Ever- 
glades, which  regulates  water  chemistry  in  the  bay.     Based  on  radiocarbon 
dating  of  mounds,   embodied  energy  in  mud  mounds  was   found  to  be  7.5  x 
1010  solar  Cal/m2. 


Energy  and  Phosphate  Landscape  of  Central  Florida 

Whereas  the  energy  systems  developing  land  forms  in  the  section  pre- 
viously considered  were  largely  without  humans,  the  world's  landscapes  are 
increasingly  dominated  by  human  management  with  huge  energies   from  fuels, 
nuclear  energy,  etc.  reforming  land.     The  study  results  of  human-nature 
interactions   in  landforms  affected  by  phosphate  strip  mining  are  given  in 
this  section.     Data  on  the  work  of  nature  in  succession  on  old  mined  land 
are  also  given  and  compared  with  human  effort  of  reclamation. 
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Embodied  Energy  in  Spoil  Mound  Form 

A  model  of  a  spoil  mound  is  shown  in  Figure  16.     The  main  energy  in- 
puts come  from  the  dragline  (controlled  by  mining  decision)   and  natural 
energies  such  as  sun,  wind,  seeds,  and  rain.     Components  of  the  spoil 
mound  are  the  ecosystem  growing  on  the  mound  and  land  area,  which  includes 
mound  form  and  small  gullies.     After  initial   input   from  the  dragline, 
mound  form  interacting  with  rain  generates  small  gullies  down  the  face  of 
the  mound  surface  through  erosion.     These  gullies  may  be  important  in 
providing  germination  sites  with  ameliorated  microclimate  for  colonizing 
vegetation.     Gullies  seem  to  accelerate  initial  development  of  ecosystem 
biomass,  which  later  spreads  outside  the  gullies  to  the  rest  of  the  mound 
surface.     The  biomass  in  turn  controls  runoff  and  thus  stabilizes  the 
system.     The  embodied  energy  of  the  mound  land  form  is  used  over  long 
periods  of  time,   interacting  with  the  succession  of  ecosystem  to  build  and 
maintain  structure. 

Embodied  energy  of  a  spoil  mound  was  calculated  as  a  by-product  of 
the  mining  process.     The  energy  signature  of  a  spoil  mound  in  the  central 
Florida  phosphate  district  is  given  in  Table  22.     The  largest  energy 
inflow  is   from  the  energy  used  by  the  dragline  in  mining  out  the  phosphate 
matrix.     Energy  input  to  the  dragline  is  30  times  greater  than  the  natural 
energy  signature,  making  it  the  dominant   feature  of  the  landscape.  The 
role  of  the  dragline  is  as  a  large  consumer — reorganizing,  recycling,  and 
even  creating  new  land.     Since  the  main  purpose  of  the  dragline  is  mining 
and  not  the  construction  of  spoil  mounds,  the  mounds  are  actually  a 
by-product  of  mining.     The  embodied  energy  used  to  mine  the  matrix  and 
thus  build  spoil  mound  structure  was  found  to  be  2.0  x  10^  solar 
Cal/mound  or  1.7  x  107  solar  Cal/m2. 
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Table  22.     Spoil  mound  case  study  energy  signature  (data  on  actual  energy 
flows  are  given  in  Appendix  II-J  and  footnotes  below). 


Energy  Source 

Actual  Energy 
Flow, 
Cal/m2-yr 

Energy  Trans- 
formation Ratio, 
solar  Cal/Cal 

Embodied  Energy 
Flow, 

106  solar  Cal/m2-yr 

Sun 

1,500,000 

1 

1.5 

Wind 

600 

74 

0.04 

Chemical  potential 
of  rain 

3.9 

3,600 

0.01 

Elevated  potential 
of  runoff 

1.5 

3,000 

0.004 

Kinetic   impact  of 
rain 

9.0 

240,000 

2.  1 

Dragline 

68.  0a 

Total  embodied  flow 

68.6 

aInputs  to  dragline  to  build  a  spoil  mound:  Inputs  to  a  35  cubic  yard 
capacity  dragline  are  given  below  from  Zeindler  (1964)  on  a  cubic  yard 
basis . 


Labor  =  2.6  x  10  3  man-hours  (at  $3.17/hr),  or  $8.2  x  10"3; 
Electricity  =  0.6  KWH; 
Maintenance  =  $0. 02; 

Amortized  capital  cost  =  $9.4  x  10~3  ($2  x  106  over  20 
years) . 

Dragline  inputs  are  about  $0.04/yd3  for  dollar  costs  and  0.6  KWH/ yd3 
for  electrical  energy  input.     Converting  these  to  solar  energy 
equivalents  and  summing  gives: 

Energy  =  ($0 . 04/yd3 ) (2  2  x  107  solar  Cal/$) 

+   (0.6  KWH/yd3)(860.5  Cal/KWH) (8000  solar  Cal/Cal) 
=  8.8  x  10*  solar  Cal/yd3  =  4.1  x  106  solar  Cal/yd3 
=  5.0  x  10b  solar  Cal/yd3. 

Using  the  average  dimensions  of  a  spoil  mound   (T.   Oxford,   pers . 
comm.)— Height  =  5  m;  width  =  20  m;  and   length  =  60  m— the  volume  of  the 
mound  can  be   found  by  assuming  prism  shape. 

Volume  =  (1/2)  height  x  width  x  length 
=  (0.5)(5  m)(20  m)(60  m) 
=  3000  m3,  or  about  4000  yd3. 
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Table  22.  (continued.) 


Total  input  to  dragline  in  building  a  spoil  mound  is  then 

Total  input  to  dragline  =  (5.0  x  loj  solar  Cal/yd3) (4000  yd3/mound) 
=  2.0  x  1010  solar  Cal/mound. 

Expressed  on  an  area  and  volume  basis  this  is 

Unit  area  input  to  dragline  =  (2.0  x  1010  solar  Cal/mound) /( 1200  m2 /mound) 
=  1.7  x  107  solar  Cal/m2. 

Unit  volume  input  to  dragline  -  (2.0  x  1010  solar  Cal/mound) /( 3000  m3 /mound) 
-  7.0  x  106  solar  Cal/m3. 

Total  overburden  and  matrix  removed  in  the  annual  operation  of  a  typical 
mine  is  22  million  cubic  yards  (Wang  et  al .  1974).  Total  energy  to  move 
this  volume  prorated  over  area  of  the  mine  is 

Dragline  energy  =   (5.0  x  106  solar  Cal/yd3) (22  x  106  yd3/yr) 
(1/400  acres  ) (acres/4047  m2) 
=  6.8  x  107  solar  Cal/m2-yr. 
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Another  aspect  of  spoil  mounds  that  affects  the  energetics  of  mining 
landscapes  is  their  surface  area.     Premining  topography  was  to  a  large  ex- 
tent flat  while  postmining  surfaces  are  diverse  with  mounds  and  depres- 
sions.    Surface  area  per  unit  of  horizontal  distance  is  greater  in  the 
mound  landscape  compared  to  the  flat  premining  landscape.     In  terms  of 
surface  area,  a  single  spoil  mound  covers  1341.6  m2  compared  to  a  cor- 
responding flat  area  of  1200  m2  (see  footnote  below).     In  other  words, 
there  is  about  10%  more  surface  area  after  mining.     The  action  of  mining 
created  new  area,  which  may  have  many  ramifications  in  terms  of  energy 
flows  that  are  input  on  an  area  basis. 

Succession  on  Overburden  Spoil  Mounds 

Development  of  Vegetation  Structure 

Results  of  measurements  of  vegetation  on  spoil  mounds  are  shown  in 
Figure  17.     A  general   increase  in  vegetation  characteristics  with  time  is 
evident  in  each  case.     Percent  cover  data  are  shown  in  Figure  17a. 
Between  1-  and  5-year-old  sites,  cover  increases  dramatically  from  about 
16%  to  about  47%.     Little  increase  was  found  from  5  to  10  years.  Cover 
increases  to  about  70%  at  the  20-year-old  sites  with  bare  area  dispersed 
amongst  grass  stems  and  clumps  rather  than  in  large  bare  patches  as  at 


Surface  area  of  spoil  mound:  Assuming  the  mound  has  the  geometry  of  a 
prism,   the  hypotenuse  of  the  right  triangle  representing  one  half  of  the 
mound  cross  section  is  h2  =  a2  +  b2  where  a  =  5  meters;  b  =  10 
meters;  and  h  =  (52  +  102)0-5  =  11.18  meters.     Surface  area  of  a 
mound  slope  is  then  A  =  (60  m)(11.18  m)  =  670.8  m2 .     Since  the  mound  has 
two  slopes  as  a  prism,   total  mound  surface  area  =  2A  =  2(670.8  m2)  = 
1341.6  m2. 
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younger  sites.    On  the  oldest  sites  cover  approaches  100%  as  indicated  by 
aerial  photographs.     Optical  density  is  shown  in  Figure  17b.     The  curve 
increases  from  0.19  to  1.02  at  San  Felasco.     This  indicates  increasing 
thickness  of  the  canopy  through  succession.     Variation  was  great  at  young 
sites  due  to  the  patchy  distribution  of  cover  but  decreased  with  site  age 
as  cover  became  more  continuous.     Tree  height  data,   shown  in  Figure  17c, 
increases  from  about  3  meters  on  the  10-year-old  site  to  a  maximum  of 
greater  than  25  meters  at  San  Felasco.     The  increase   is  nearly  linear  with 
no  indication  of  leveling  at  the  oldest  sites. 

Development  of  Litter  and  Soil 

Figure  18  presents  data  on  litter,  soil  organic  matter,  and  ants  as 
an  index  of  soil-animal  activity.     The  general  trend  again  is  an  increase 
over  time  in  each  case.     Litter  weight,   shown  in  Figure  18a,  develops 
quickly  during  the  first  10  years,  slows  between  years  10  and  20,  then 
increases  in  later  stages  to  a  maximum  of  almost  1600  g/m2.  Litter 
depth  or  thickness  data  are  given  in  Figure  18b.     Litter  was  too  sparse  to 
measure  depth  at  the  1-year-old  site.     Litter  depth  drops  slightly  between 
5  and  10  years  then  increases  gradually  to  the  50-year  level.     A  large 
increase  was  found  between  the  50-year-old  sites  and  San  Felasco,  propor- 
tionally greater  than  the  increase   in  litter  weight.     This   indicates  that 
the  litter  at  San  Felasco  has  relatively  more  air  space  than  at  the 
50-year-old  site.       Figure  18c  shows   the  development  of  soil  organic 
matter.     After  an  initial  drop  between  1  and  5  years,  soil  organic  matter 
increases   to  a  maximum  of  over  8%  at  the  50-year-old  sites   then  drops  at 
San  Felasco.     Ant  nest  data  for  the  first  10  years  are  shown  in  Figure 
18d.     The  increase   in  ant  nests  over  time  is  dramatic,   from  2  per  10  m2 
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at  1  year  to  over  8  per  10  m2  at  the  10-year-old  site.  Also  a  general 
increase  in  ant  nest  size  with  increasing  site  age  was  noted  up  to  the 
50-year-old  sites  where  few  ants  were  seen. 

Contrasts  in  soil  profile  development  are  shown  in  Figure  19.  The 
1-year-old  soil  has  a  hardened  surface  with  little  differentiation  down 
the  profile.     The  50-year-old  soil  shows  considerable  development  relative 
to  the  1-year-old  soil.     Figure  19  shows  several  strata  evident  including 
litter  already  mentioned.     Roots  and  organic  matter  are  interspersed 
throughout  the  upper  layers. 

Forest  Structure  and  Composition  on 
50-Year-01d  Sites  and  San  Felasco 

Basic  data  of  stem  density  and  basal  area  for  tree  species  plus 
standing  dead  at  the  two  50-year-old  sites  are  compared  in  Table  23.  The 
same  species  are  found  at  both  sites  with  moderate  changes  in  characteris- 
tics between  sites.     Importance  values  are  calculated  in  Table  24.  Both 
stands  are  dominated  by  live  oak,  Quercus  virginiana,  with  subdominates 
switching  importance  between  sites. 

Tables  25  and  26  give  data  for  canopy  trees   (greater  than  10  centi- 
meters dbh)  at  San  Felasco  and  the  phosphate  spoil  mounds,  respectively. 
Basal  area  and  stem  density  are  combined  as  importance  values   for  each 
species.     At  San  Felasco,   17  species  (counting  standing  dead)  were  found 
with  a  total  basal  area  of  35.4  m2/ha  and  a  total  stem  density  of  602 
trees/ha.     Six  species  (counting  standing  dead)  were   found  on  the  two 
50-year-old  phosphate  spoil  mounds  with  a  total  basal  area  of  28.6  m2/ha 
and  a  total  stem  density  of  368  trees/ha. 
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Table  23.     Comparison  of  tree  data  at  the  50-year  sites. 


Payne  Creek  Christ ina 


Stem  Basal  Stem  Basal 

Density,       Area,  Density,  Area, 

trees/ha      m^/ha  trees/ha  m^/ha 


Quercus  virginiana  Mill. 

122 

11.4 

135 

21.9 

Quercus  nigra  L. 

52 

4.0 

74 

6.5 

Liquidambar  styraciflua  L. 

93 

4.0 

48 

2.3 

Pinus  elliottii  Engelm. 

70 

2.0 

9 

0.4 

Prunus  serotina  Ehrh. 

7 

0.6 

61 

2.2 

Standing  dead 

19 

0.4 

48 

2.4 

Total 

363 

22.4 

375 

35.7 
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Table  25.     San  Felasco  Hammock  tree  data.     Numbers  in  parentheses  are  per- 
centages . 


Stem 
Density, 
trees/ha 

Basal 
Area, 
m2/ha 

Importance 
Value 

Magnolia  grandiflora  L. 

100 

(16. 

6) 

9 . 

8 

(27. 

7) 

22 .  15 

Liquidambar  styraciflua  L. 

137 

(22. 

8) 

7. 

0 

(19. 

8 ) 

21 . 30 

Tilia  floridana  Small 

63 

(10. 

5) 

4  t 

1 

(11. 

6) 

11.05 

Quercus  laurifolia  Michx. 

44 

(7. 

3) 

4. 

2 

(11. 

9) 

9.  60 

Ostrya  virginiana  (Mill.) 

K.  Koch 

78 

(13. 

0) 

1. 

7 

(4. 

8) 

8.90 

Quercus  marilandica  Muench  H. 

19 

(3. 

1) 

4. 

5 

(12, 

.7) 

7.90 

Carpinus  caroliniana  Walt. 

70 

(11. 

6) 

1 , 

2 

(3. 

4) 

7 .  5 

Magnolia  virginiana  L. 

15 

(2, 

.5) 

1 , 

0 

(2. 

8) 

2.  65 

Carya  glabra  (Mill.)  Sweet 

19 

(3. 

,1) 

0 . 

6 

(1, 

,  7) 

2 . 40 

Cornus  florida  L. 

19 

(3. 

1) 

0. 

4 

(1, 

1 ) 

2  10 

Nyssa  sylvatica  var .  biflora 

(Walt.)  Sarg. 

11 

(1. 

.8) 

0. 

3 

(0. 

8) 

1.30 

Pinus  glabra  Walt. 

4 

(0. 

7) 

0. 

3 

(0. 

8) 

0.75 

Quercus  nigra  L. 

7 

(1, 

,1) 

0. 

1 

(0. 

3) 

0.  70 

Ilex  opaca  Ait. 

4 

(0. 

7) 

0. 

1 

(0. 

3) 

0.50 

Robinia  pseudoacacia 

4 

(0. 

,7) 

>0. 

1 

(0. 

.1) 

0.40 

Acer  rubrum  L. 

4 

(0. 

7) 

>0. 

1 

(0. 

1) 

0.40 

Standing  dead 

4 

(0. 

7) 

>0. 

1_ 

(0, 

1) 

0.40 

Total 

602 

35. 

4 
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Table  26.     Combined  phosphate  spoil  mound  forest  tree  data.     Numbers  in 
parentheses  are  percentages. 


Stem 

Basal 

Density, 

Area, 

Import anc  e 

trees/ha 

m2 

Val  ue 

Quercus  virginiana  Mill. 

128 

(34.8) 

16.2  (56.8) 

45.80 

Quercus  nigra  L. 

62 

(16.8) 

5.1  (17.9) 

17.35 

Liquidambar  styraciflua  L. 

72 

(19.6) 

3.2  (11.2) 

15.40 

Pinus  elliottii  Engelm. 

42 

(11.4) 

1.3  (4.6) 

8.00 

Prunus  serotina  Ehrh. 

32 

(8.7) 

1.3  (4.6) 

6.65 

Standing  dead 

32 

(8.7) 

1 .4  (4.9) 

6.80 

Total 

368 

28.5 
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Standing  dead  trees  have  low  importance  values  but  are  a  constant 
feature  of  the  forests  as  indicated  by  the  frequency  data.     They  may  play 
important  functional  roles  in  forest  ecology  in  terms  of  wildlife  habitat 
and  nutrient  conservation  (Egler  1975;  Bormann  and  Likens  1979).  Impor- 
tance values  of  dead  trees  at  the  50-year-old  sites  was  17  times  greater 
than  at  San  Felasco,  perhaps  reflecting  that  the  tree  community  of  the 
phosphate  forest  has  gone  through  more  recent  transitions  and  is  not  yet 
at  steady  state. 

Simulation  of  Succession 

The  succession  model  with  data  used  for  evaluation  of  coefficients  is 
shown  in  Figure  20.     Storage  units  are  g  dry  weight/m2  and  flow  units 
are  g  dry  weight/m2 • yr .     References  and  derivations  are  given  in 
Appendix  III  and  the  program  in  included  as  Appendix  IV.     The  standard 
simulation  run  of  succession  is  given  in  Figure  21a.     Grasses  peak  first, 
after  about  5  years;  pines  dominate  next  from  about  10  to  40  years;  climax 
hardwoods  supplant  pines  and  reach  steady  state  at  about  60  years;  and 
litter  gradually  builds  to  a  steady  state  after  about  50  years.  This 
general  pattern  of  vegetation  change  fits  the  field  observations  summar- 
ized in  the  discussion. 

Table  27  provides  validation  of  the  model  by  comparing  predicted  lit- 
ter biomass,  run  at  different  time   integration  intervals,  with  litter  bio- 
mass  from  field  sampling  of  phosphate  mine  ecosystems  (see  Figure  18b). 
With  time  steps  of  0.25  years  and  0.1  years,  only  the  predicted  5-  and 
50-year  litter  values  fall  within  the  standard  error  limits  of  observed 
litter.     With  a  time  step  of  0.025  years,   the  predicted  20-year  litter 
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Figure  21.     Comparison  of  the  effects  of  increasing  the  seeding  rate  on  the 

rate  of  succession,     a.  standard  run,  b.  result  of  increasing  the 
seeding  rate  1000  times  after  the  litter  germination  thresholds 
are  reached.     Grass  vegetation  biomass,  Q^,  is  scaled  by  4  times 
to  excentuate  its  growth  curve. 
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Comparison  of  litter  values  from  simulation  runs  with  varying  time 
steps  (DT),  with  litter  data  from  field  sampling  for  validation. 


r>  Measured 
yr  Q4(DT=0.25)      Q4(DT=0.1)       Q4(DT=0.025)      Q4(DT=0.01)  Litter 


1 

4.0 

12.4 

36.6 

61.8 

54.0 

+ 

11. 

8 

5 

309.0 

317.5 

340.0 

375.2 

377.  2 

+ 

61, 

.6 

10 

341.0 

342.3 

369.7 

342.8 

622.8 

71. 

5 

20 

524.5 

557.2 

604.3 

668.9 

688.  5 

96. 

7 

50 

1196.8 

1223.13 

1257.4 

1300.5 

1250.4 

+ 

73. 

1 
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also  fits  the  data.  Best  fit  of  predicted  values  with  those  observed  was 
found  with  a  time  interval  of  0.01  years,  for  which  all  predictions  except 
the  10-year  value  fell  within  acceptable  variation.  Due  to  the  magnitude 
of  the  scaling  factors  and  length  of  simulation  runs  (70  years),  however, 
these  changes  in  precision  did  not  greatly  effect  the  form  of  the 
simulation  graphs.  Therefore,  simulation  graphs  were  plotted  with  a  time 
step  of  0.25  years  since  this  ran  much  faster. 

Two  principal  types  of  simulations  were  performed  with  the  model: 
increasing  seeds  as  a  pulse  after  germination  thresholds  were  reached  and 
adding  litter.     The  performance  of  these  experiments  on  the  model  are  sum- 
marzied  in  Table  28.     Performance  was  judged  by  the  age  at  which  hardwood 
steady  state  was  reached.     All  simulations  with  seed  additions  had  the 
same  pattern  of  growth  and  transition  of  vegetation  as  the  standard  or 
control  run  shown  in  Figure  21a.     Succession  was  accelerated  from  5%  to 
25%  by  this  type  of  simulation.    As  an  example  the  run  with  all  seeds 
being  pulsed  1000  times  is  shown  in  Figure  21b.     Relative  to  the  control 
run,  the  first  two  stages  in  succession  are  abbreviated  even  though  they 
still  reach  comparable  peak  biomass  values  as   in  the  control  succession. 
Climax  is  reached  in  year  44,  which  is  12  years  faster  than  in  the 
control . 

The  effects  of  adding  litter  with  and  without  light  threshold  on 
hardwood  germination  are  shown  in  Figure  22a  and  22b,  respectively,  and 
also  in  Table  28.     With  the  light  control,  no  addition  of  litter  caused 
significant  change  in  succession,   since  hardwood  germination  was  being 
limited  by  light,  not  litter.     Tests  run  on  the  model  without  the  light 
control  showed  that  the  optimal  litter  addition  was  300  g/m2  on  year  5 
of  succession.     This  strategy  allowed  the   fast  growing  grass  stage  to 
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Table  28.     Comparison  of  performance  of  simulation  experiments  in  acceler- 
ating succession. 


Simulation  Age  of  Percent 

Description  Climax,   yr  Acceleration 


Standard  run  59 

Pulse  grass  seed  by  1000  56  5 

Pulse  pine  seed  by  1000  54  9 

Pulse  hardwood  seed  by  1000  52  12 

Pulse  pine  and  hardwood 

seeds  by  1000  47  20 

Pulse  all  seeds  by 

1000  44  25 

100  51  14 

10  56  5 

Add  litter  w/light  control  59  0 

Add  litter  w/ o  light  control  41  31 

Pulse  all  seeds  by  1000  and 

add  litter  33  44 
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Figure  22.     Comparison  of  the  effects  of  the  light  threshold  on  litter 

additions  to  the  simulation  model  of  succession,  a.  addition 
of  litter  with  the  light  control  threshold,  b.  addition  of 
litter  without  the  light  control  threshold.  Grass  vegetation 
biomass,  Q^,  is  scaled  by  4  times  to  excentuate  its  growth  curve. 


develop  its  peak  litter  biomass  to  which  the  minimum  amount  of  litter  was 
added  to  initiate  germination  of  hardwoods.     This  experiment  shortened  the 
pine  stage  and  decreased  its  peak  biomass,  perhaps  since  its  role  in  suc- 
cession was  taken  up  in  part  by  the  litter  addition.     Climax  is  reached  in 
year  41,  which  is  3  years  faster  than  in  the  seeding  experiments  and  18 
years  faster  than  the  control  simulation. 

A  final  simulation  combined  the  addition  of  seeds  and  litter  without 
the  light  control  on  hardwood  germination  to  view  the  maximum  amount  that 
succession  may  be  able  to  be  accelerated.     This  run  is  shown  in  Figure 
23b.     The  general  pattern  of  succession  is  similar  to  that  found  with  the 
addition  of  litter  (Figure  22b).     However,  climax  is  reached  by  year  33, 
which  is  44%  faster  than  in  the  standard  or  control. 

Energy  Analysis  of  Reclamation  Alternatives 

Many  unreclaimed  spoil  mounds  of  various  ages  are  scattered  through- 
out the  central  Florida  phosphate  district  left  from  past  mining.  Recla- 
mation of  these  and  recently  produced  mounds  often  consists  of  leveling 
the  spoil  mounds  and  planting  pasture  grasses.     However,   in  many  cases 
succession  has   filled  the  role  of  reclamation  through  the  development  of 
ecosystems  on  the  abandoned  mines.     An  important  issue  for  state  and  local 
governments  is  to  decide  whether  these  old  mined  lands  need  to  be  re- 
claimed.    One  way  this  problem  can  be  approached  is  by  comparing  embodied 
energy  in  the  two  reclamation  alternatives. 

A  phosphate  mining  model   is  shown  in  Figure  24.     This  model  does  not 
describe  the  entire  phosphate  production  industry  but  emphasizes  land  rec- 
lamation alternatives,   including  putting  land  into  pasture  through 
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Figure  23.     Comparison  of  the  effects  of  increasing  the  seeding  rate  and  addi- 
tion of  litter  on  the  rate  of  succession,  a.  stand  run;  b.  result 
of  increasing  seeding  rate  1000  times  and  adding  litter.  Grass 
vegetation  biomass,  Q^>  is  scaled  by  4  times  to  excentuate  its 
growth  curve. 
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conventional  reclamation  versus  natural  succession  restoring  the  ecosys- 
tems.    The  extraction  of  the  phosphate  matrix  is  shown  as  being  exported 
for  further  processing.     The  model  is  based  on  the  annual  operations  of  a 
representive  mine  with  400  acres  of  active  mining.     Embodied  energy  flows 
are  shown  as  1012  solar  Cal/mine-yr.     The  results  of  the  calcula- 
tions are  summarized  in  Table  29. 

The  embodied  energy  equations  describing  the  two  reclamation  alterna- 
tives are  set  up  as  an  equality  in  Figure  25.     The  age  of  spoil  mound  that 
satisfies  the  equality  is  then  solved  for  by  rearranging  the  equations  and 
substituting  the  calculated  values   from  Table  29.     This  solution  repre- 
sents the  threshold  age  of  an  unreclaimed  mound  that  has  more  embodied 
energy  from  succession  than  could  be  generated  by  mechanical  reclamation. 
This  threshold  was  found  to  be  at  age  16  years.     Conventional  reclamation 
of  mined  lands  of  this  age  and  younger  could  yield  an  embodied  energy  sav- 
ings, but  on  older  lands  it  would  yield  less  embodied  energy  than  it  con- 
sumes.    If  the  end  land  use  is  taken  to  be  improved  pasture  the  cost  of 
leveling  the  spoil  mounds  is  $750/acre  (Wang  et  al .   1974).    Without  mining 
the  cost  of  clearing  and  preparing  the  land   for  pasture  is  $136/acre  (see 
footnote  below).     Therefore  more  than  $600/acre  in  additional  costs  are 
required  to  develop  pasture  on  mined  land  versus  unmined  land.     This  dif- 
ference makes  the  payback  period  longer  for  cattle  operation  on  reclaimed 
mined  lands.     However,    improved  pasture  may  be  the  best   landuse   in  the 
long  run  for  reclaimed  land  from  the  raining  industry's  perspective  since 
it  generates  economic  return.     Natural  succession  generates  no  direct 


Cost  of  clearing  and  preparing  land  for  improved  pasture:  Bulldozer 
costs,   $542/day  (Godfrey  1976)  divided  by  4  acres  cleared/day  (Ruess 
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Table  29.     Embodied  energy  flows   in  Figure  24. 


Flow 

Embodied  Energy 
10*2  solar  Cal/mine'yr 

Footnote 

Natural   input  (sun,  wind,  rain) 

3.4 

a 

Lost  productivity 

2.2 

b 

Energy  in  reclamation 

6.6 

c 

Inputs  to  dragline 

110.0 

d 

Pasture  subsidy 

0.  5 

e 

'Natural  energy  input:     The  natural  energy  portion  of  the  energy  signa- 
ture for  a  spoil  mound   is  2.1  x  106  solar  Cal/m2-yr  (Table  22). 
Expressed  over  the  area  of  the  mine  the  energy  equivalent  is 

Energy  =  (2.1  x  106  solar  Cal/m2 • yr ) (400  acres) (4047  m2/acre) 

=  3.4  x  1012  solar  Cal/yr. 


Energy  in  lost  productivity:     As  a  consequence  of  mining,  gross  primary 
productivity  is  lost  over  the  length  of  time  needed  to  return  to  past 
levels.     Using  productivity  of  pineland  from  Costanza  (1975)  and  approxi- 
mating total  loss  as  one-half  of  this  rate  for  20  years  of  succession,  the 
energy  equivalent  is 

Energy  =  (5.5  x  106  coal  Cal/acre -yr) (2000  solar  Cal/coal  Cal) 
(400  acres)(0.5) 

=  2.2  x  1012  solar  Cal/yr. 

cEnergy  in  reclamation:     The  costs  of  phosphate  mine  reclamation  are 
$750/acre  (Wang  et  al.   1974).     The  energy  equivalent  expressed  over  the 
area  of  a  typical  mine  is 

Energy  =   ($750/acre) (2 . 2  x  107  solar  Cal/$)(400  acres/yr) 

=  6.6  x  1012  solar  Cal/yr. 

dDragline  inputs:     Dragline  inputs  were  calculated   in  the  footnote  to 
Table  22  at  5  x  106  solar  Cal/yd3.     Total  overburden  and  matrix  re- 
moved in  the  annual  operation  of  a  typical  mine  is  22  million  yd3  (Wang 
et  al.     1974).     Total  energy  equivalent  is 

Energy  =   (22  x  1 06  yd3/yr)(5  x  106  solar  Cal/yd3) 

=  1.1  x  1014  solar  Cal/yr. 


Table  29. 


( cont  inued . ) 


ePasture  subsidy:     Costanza  (1975)  gives  the  fossil  fuel  subsidy  to 
improved  pasture  at  6.29  x  105  coal  Cal/acre-yr.     Expressed  for  the 
mine  area  the  energy  signature  is 

Energy  =  (6.29  x  105  coal  Cal/  acre *yr) (2000  solar  Cal/coal  Cal) 
(400  acres) 


=  0.5  x  1012  solar  Cal/yr. 
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value  but  its  indirect  values  of  restoring  fertility  and  as  wildlife  habi- 
tat may  be  significant. 

Phosphate  Fertilizer  Industry  Model 

A  model  of  the  Florida  phosphate  fertilizer  industry  is  shown  in  Fig- 
ure 26.     The  overall   industry  is  depicted  as  a  main  circuit  of  phosphate 
flowing  from  the  landscape  to  the  fertilizer  market.     The  process   is  init- 
iated and  maintained  by  a  switch,  which  senses  the  amount  of  phosphate  in 
the  matrix.     Mining  will  only  take  place  if  sufficient  deposits  are  antic- 
ipated.    The  industry  is  shown  with  two  main  parts,  the  mine  itself  and 
the  fertilizer  plant.     Economic  exchange  is  with  the  fertilizer  market  in 
the  outside  economy.     Inputs  to  the  industry  include  water,  matrix, 
capital   investment,   and  goods,   services,   and  fuels.     Water  is  mostly  re- 
cycled.    The  mine  includes  strip  raining  by  draglines,  which  remove  over- 
burden spoil  and  expose  the  matrix,  water   flows,   which  are  used  to  trans- 
port the  matrix  and  wastes  as  a  slurry  in  pipes,  and  the  beneficiation 
plant  where  matrix  is  separated  into  sand,   clay,   and  phosphate,   which  is 
sent  to  the  fertilizer  plant.     There  the  phosphate  is  mixed  with  sulfuric 
acid,   upgrading  the  phosphate  rock  to  produce  commercial  fertilizer. 
Gypsum  is  a  by-product  waste  of  this  process.     Reclamation  is  shown  as  an 
"upstream"  circuit  restoring  spoil,   sand  tailings,   and  slime  clays.  This 
includes  land  manipulation,  such  as  grading  and  diking,  plus  some  revege- 
tation.     Gypsum  is  not  reclaimed,  but   is   left   stacked  next  to  the  fertil- 
izer plant. 

An  aggregated  model  of  the  phosphate  fertilizer  industry  is  shown  in 
Figure  27.     This  model   is  evaluated   for  cost  of  mining  and  processing 


phosphatic  matrix  into  phosphoric  acid   fertilizer,  on  the  basis  of  a 
single  typical   industry  unit.     Total  costs  are  4.6  x  1015  solar  Cal/yr 
for  an  annual  fertilizer  output  of  about  1.2  x  106  tons.     In  the  annual 
production  of  fertilizer,  2.4  x  107  tons  of  waste  land  material  (about 
20  tons  of  waste/ton  fertilizer)   are  created  (see   footnote  below). 
Reclamation  of  these  waste  lands  requires  1.4  x  1013  solar  Cal.  This 
model   illustrates  the  relative  magnitude  of  reclamation  in  terms  of  the 
total  budget.     Reclamation  is  0.3%  of  the  total,  showing  that  its  overall 
cost  is  very  small.     Calculations  of  these  numbers  are  given  in  Appendix 
V. 

Essentially  all  of  the  embodied  flow  of  fertilizer  to  market  is  net 
energy  (55.8  -  4.6),  which  has  a  yield  ratio  of  12  solar  Cal  output  per 
solar  Cal  of  input.     These  large  values  are  the  result  of  mining  a  rich, 
nonrenewable  storage. 


Waste  volumes  generated  in  the  production  of  fertilizer.     Spoil : 
13,000,000  yd3  (EPA  1978)/5  x  106  ton  of  rock  (Greene  1978)  =  276~~y73/ton 
of  rock;   (2.6  yd3/ton  of  rock)(3.31  ton  of  rock/ ton  of  fertilizer) (Phil- 
lips 1975)  =  8.6  yd  /ton  of  fertilizer.     Assume  density  of  spoil  as  1.4 
g/crn3.     [(8.6  yd3/ton) (0. 7646  m3/yd3)(106  cm3/m3)(1.4  g/ cm3) (kg/ 1000  g)]/ 
(10    kg/ ton)  =9.2  tons.     Clay:  One  ton  of  clay  is  removed   for  each  ton 
of  rock  (EPA  1978)    (1  ton  clay/1  ton  rock) (3.31  ton  rock/ ton  of  fertil- 
izer) =  3.31  tons.     Sand:   Same  as  clay  (EPA  1978);   3.31  tons.     Gypsum:  4.5 
tons  (Bowers  1968).     Total  wastes  =  20.32  tons/ton  of  fertilizer.  Total 
production  of  fertilizer  is  1.2  x  106  tons/yr,  making  production  of 
wastes  =  (1.2  x  10b  tons/yr) (20 . 32  tons/ton)  =  2.4  x  107  tons/yr. 


CHAPTER  4 
DISCUSSION 


Energy  Causality  of  Landscape 

Energy  has  been  suggested  as  the  basis  for  causal  action  in  all 
levels  and  kinds  of  real  systems.     Theories,  models,  and  calculations  have 
been  given  here  for  the  energy  causality  of  landscape.     A  theory  of  the 
spatial  expression  of  energy  sources  and  the  complementary  shapes  they 
cause  to  develop  in  landscape  components  is  central  in  these  explanations. 
The  theory  as  given  here  states  that  the  pattern  of  incoming  energies 
interacting  with  a  surface  builds  forms  on  the  surface,  which  amplify  the 
spatial  program  of  energy.     This  amplification  is  expressed  in  the  conver- 
gence and  divergence  of  energies  as  they  flow  across  a  landscape.  For 
example,   the  form  of  a  valley  causes  runoff  on  its  slopes   to  converge  into 
a  central  river  channel.    Without  the  valley  form  the  power  in  runoff 
would  be  less   focused  or  concentrated  and  thus  the   form  has  amplified  the 
energy  flow. 

Because  landscape  itself  is  a  map  of  accumulated  effect  of  its  his- 
torical energy  inputs  and  transactions,  the  patterns  can  reveal  something 
about  how  energies  are  used.  One  hypothesis  is  that  the  proportion  of  a 
landform  type  is  directly  related  to  the  proportion  of  its  matching  spa- 
tial energy  type  in  the  energy  signature  of  the  landscape.  This  hypoth- 
esis is  tested  by  calculating  energy  signatures  for  three  landscapes  dom- 
inated by  single  landform  types  and  one  landscape  with  a  complex  pattern 
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of  overlapping  landform  types.     If  the  energy  input  algorithm  and  the 
energy  quality  factors  are  correct,   for  the  first  three  landscapes  the 
dominant  spatial  energy  type  should  match  with  that  indicated  by  the  domi- 
nant landform  type.     In  the  fourth  case,  no  single  spatial  type  of  energy 
should  predominate. 

The  landscapes,  all  from  Florida,  have  been  studied  with  construction 
of  maps  and  energy  signatures.     The  Green  Swamp  area  in  central  Florida  is 
dominated  by  cypress  domes,  which  are  circular  to  irregularly  shaped 
depression  wetlands  scattered  in  pine  flatwoods.     Brown  (1978)   and  Odum 
(1978)  discuss  the  role  of  these  wetlands  in  the  landscape.     In  this  exam- 
ple, the  cypress  domes  are  predicted  as  center  shapes  generated  by  conver- 
gence of  sheet  energies  from  the  background  of  pine  flatwoods.     A  sheet 
energy  should  then  dominate  this  energy  signature.     An  energy  signature 
for  the  Green  Swamp  is  given  by  McCuller  (1975).     The  Apalachicola  River 
basin  is  located  in  the  Florida  panhandle  and  extends  north  into  Georgia. 
The  dominant  feature  of  this  landscape  is  the  floodplain  forest  along  the 
river  and  is  reflective  of  line  energy,  which  dominates  this  energy  signa- 
ture.    An  energy  signature  for  the  Apalachicola  region  is  given  by  Boynton 
(1975).     The  Florida  coral  reef  is  a  barrier  reef  extending  along  the 
Atlantic  side  of  the  Florida  Keys.     This  is  predicted  as  a  zone  shape  and 
the  energy  signature  is  shown  in  Table  8  (see  page  49).     The  dominant 
energy  in  this  energy  signature  should  be  frontal  type.     Hendry  County,  in 
southwest  Florida,   is  the  final   landscape  considered.     It   is  dominated  by 
a  variety  of  wetlands:     cypress  strands  and  plains,  extensive  marshes  bor- 
dering lake  Okeechobee,   and  scattered  seasonal  wet  prairies  and  marshes. 
This  is  used  as  an  example  of  a  complex  landscape  with  several  shape  types 
with  no  one  type  dominant.     Thus,   no  dominance  should  occur   in  the  energy 
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signature  by  specific  spatial  type.  An  energy  signature  for  Hendry  County 
is  given  by  DeBellevue  (1975). 

A  summary  of  the  energy  signatures  is  given  in  Table  30  by  their 
spatial  type.     The  predictions  hold  for  the  Apalachicola  Basin  and  the 
Florida  coral  reef.     Line  dominance  is  80  times  for  the  floodplain  land- 
scape and  frontal  dominance  is  more  than  1000  times  for  the  coral  reef. 
For  the  Green  Swamp,  a  sheet  energy  is  the  largest  energy  but  it  is  only 
slightly  larger  than  the  maximum  frontal  energy.     For  Hendry  County  maxi- 
mum sheet  energy  was  about  2  times  as  large  as  the  maximum  frontal  energy, 
which  may  be  similar  enough  in  magnitude  to  mix  form  types   into  a  compos- 
ite. 

Most  landscape  are  combinations  of  the  basic  form  types  blended 
together  with  no  single  type  dominant  as  with  Hendry  County.     In  these 
cases  the  ratio  of  form  types  may  not  be  directly  proportional  to  the 
ratio  of  spatial  energies  but  may  be  more  a  function  of  the  way  the  ener- 
gies interact.     Unraveling  composite  patterns  is  the  key  to  understanding 
landscape.     Comparative  study  of  different  landscapes  and  their  energy 
signatures  may  reveal  interaction  functions  of  different  energies. 

Energy  and  the  Shape  of  Landforms 

Land  form  has  been  defined  as  an  energy  storage  of  landscape  in  terms 
of  volume,  either  in  depressions  or  elevations.     This  was  quantified  by 
knowing  the  relative  dimensions  of  the  volume  and  by  assuming  some  geomet- 
ric shape.     Why  are  landform  shapes  different?     It   is  proposed  here  that 
energy  is  used  more  effectively  when  shape  matches  the  spatial  expression 
of  energy  input.     Energy  builds  the  structures  and  thus  controls  the  shape 
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Table  30.     Maximum  embodied  energy  within  energy  signatures  of  four  Florida 
landscapes  according  to  spatial  expression.     Values  are  ex- 
pressed x  10^  solar  Cal/m    yr . 


Green 
Swamp 

Apalach  i- 
cola  Basin 

Florida 
Coral  Reef 

Hendry 
County 

Maximum 

sheet  energy3 

1.4 

1.4 

1.4 

1.4 

Maximum 

frontal  energy 

0.9b 

3.7C 

1800d 

0.6e 

Max  imum 

line  energy 

310f 

aSunlight . 
bWind. 

cPhysical  runoff. 
dWaves. 

eChemical  runoff, 
f River  head. 


Ill 

of  the  structures.     Special  shapes  may  aid  in  capturing,  absorbing,  and 
possibly  using  energy  inputs.     These  shapes  are  believed  to  be  adaptive  to 
specific  spatial  energy  distributions.     Evidence  in  the  preceding  section 
in  part  supports  this   idea,  when  one  spatial  energy  type  dominated  an 
energy  signature  the  corresponding  shape-type  dominated  the  landscape.  If 
a  certain  shape  allows  more  control,   as  measured   in  Figure  1   (see  page  2) 
by  the  feedback  pathway  into  the  production  function,  it  can  be  "selected 
for"  in  competition  with  landform  alternatives  of  other  shapes. 

Shapes  of  landforms  can  aid  in  concentrating  energies  such  as  flows 
of  water,  wind  in  desert,   and  waves  and  tides  on  coasts.     Sheets  can  be 
concentrated  to  form  points,  and  fronts  can  be  concentrated  to  form  lines. 
Once  concentrated,   the  new  energies   interact  and  converge  other  less  con- 
centrated energies  from  their  surroundings.     Energy  lines  cause  a  linear 
convergence  and  energy  points  cause  a  radial  convergence,  which  helps 
maintain  the  structure  they  build.     If  the  concentrated  point  or  line  does 
not  find  enough  energy  with  which  to  match  or  interact,   it  will  disperse. 
For  example,  consider  the  upper  end  of  a  watershed  catchment  where  runoff 
moving  as  a  front  coverges  to  form  a  concentrated  line   in  the  form  of 
river  flow.     The  line  collects  more  energy  as  it  proceeds  downslope  until 
it  reaches  the  coastal  plain  where  the  energy  gradient  decreases  and  the 
line  must  disperse  back  into  a  front  as  the  river  forms  a  delta  or  wide 
anastomosing  floodplain.     Concentrated  sources  can  be  insulated   for  some 
purposes,  but  often  at  great  cost.     An  example  is  a  channelized  river. 

A  classification  of  landscape  types  is  proposed  based  on  the  spatial 
expressions  of  energy  (Figure  28).  Implied  is  that  special  forms  develop 
in  response  the  particular  spatial  energy  types.  Based  on  study  of  maps, 
six  basic  categories  of  form  are  proposed:  background,  center,  zone, 
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ZONE 


STRIP 

Figure  28.     Patterns  of  structure  postulated  to  result  from  energy  inflow 
patterns.     Names  of  the  dominant  landscape  form  type  are  given 
to  the  left  (see  Table  31)  . 
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linear  strip,   transverse  string,   and  streamlined  island.     These  are  common 
shapes  of  land  surface  recognizable  on  vegetation  and  land  use  maps. 
Different  forms  are  distinguished  by  their  shape,  orientation  and  to  some 
extent  relative  size  (Table  31).     This  has  similarities  to  a  shape 
classification  of  Florida  wetlands  (Wharton  et  al .    1976),  which  recognized 
linear,  broad  basin,  and  circular  (including  irregular  types)  shapes  in 
relation  to  water  flows.     Odum  (1978a)  presented  a  related  series  of 
wetlands  whose  form  changes  along  a  gradient  of  increasing  water  flows. 

Several  hypotheses  relate  landform  to  sources  as  follows:  recurring 
patterns  of  form  arise  under  similar  external  energy  conditions.  The 
forms  are  organizational  arrangements  for  energy  use.     Spacing  may  be 
directly  related  to  the  absorption  and  dissipation  of  energy.     The  spacing 
of  individual  energy-absorbing  forms  may  be  important  in  organizing  the 
larger  patterns.     If  spacing  is  too  close,  competition  for  energy  may 
occur,  and  if  spacing  is  too  great,  some  energy  may  pass  without  use. 

Landscape  shape  also  changes  along  hierarchies  of  converging  energy. 
Central  place  type  hierarchies  are  well  known  and  represent  the  concentra- 
tion of  sheet  energy  over  broad  areas  into  center  structures  that  perform 
services  for  the  surrounding  region. 

Lateral  hierarchy  develops  where  energies  concentrate  in  relation  to 
frontal  energies.     The  sequence  in  this  hierarchy  is  suggested  in  Figure  4 
progressing,  with  increasing  energy,   from  zones  to  strips.     An  example  is 
a  watershed  where  the  slopes  capture  and  laterally  focus  water  to  a  strip, 
which  is  the  river  channel.     Strings  and  islands  fall  within  this  ideal 
hierarchy  under  certain  intermediate  conditions  of  energy.     It   is  hypothe- 
sized that  at  low  flows,   landforms  are  arranged  as  strings  that  maximize 
the  area  exposed  to  the  flows.     As  the  flow  increases,   the  strings  rotate 
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90°  to  minimize  the  area  directly  exposed  to  the  flow.     Examples  are 
desert  dunes  and  some  peat  mounds   in  wetlandscapes . 

The  storages  in  central  place  and  "lateral  place"  hierarchies  may 
have  the  same  embodied  energy  distribution  but  the  shapes  of  the  struc- 
tures differ  because  of  the  generic  difference  of  source  energy  input, 
either  parallel  or  perpendicular  to  the  surface. 

An  advantage  of  the  shape  classification  is  that  it  may  allow  a  view 
of  a  class  of  emergent  properties  not  recognized  before.     Ecologists  have 
studied  in  detail  the  ecosystem  adaptations  for  assimilating  sunlight, 
such  as  chlorophyll  content,  assimilation  number,  photic  zone  thickness, 
and  leaf  area,  some  of  which  are  recognized  as  emergent  or  summary  proper- 
ties.    However,   these  apply  only  to  one  of  the  energy  sources  utilized  by 
ecoystems  (sunlight).     Energy  signatures  are  composed  of  other  energy 
sources  (auxiliary  energies  or  energy  subsidies)   that  have  been  shown  in 
many  cases  to  be  of  greater  magnitude  than  direct  sunlight.     Landform  may 
be  a  synoptic  ecosystem  adaptation  for  assimilating  auxiliary  energies. 
The  recognition  of  different  landform  shapes  is  a  step  toward  investigat- 
ing the  ways  landscapes  capture  and  assimilate  energy  sources.  Probably 
few  land  systems  are  developed  with  sunlight  as  the  dominant  energy 
source . 

Table  32  compares  shape  types  in  four  different  kinds  of  landscape 
disciplines.     The  classification  of  landscape  components  is   from  the  pres- 
ent study.     The  mental   image  classes  come  from  Lynch  (I960),  urban  classes 
are  a  mixture  of  city  subsystems  (see  Clay  [1973]    for  more  discussion), 
and  the  components  of  art  come  in  part  from  Klee  (1953).     This  convergence 
of  thought  indicates  a  general  principle  of  the  categories  of  shapes  that 
are  adaptive  and  develop  in  many  different  kinds  of  energy  systems.  In 
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Table  32.     A  comparison  of  shape  classifications  in  ecology,  environmental  per- 
ception, urban  geography,   and  art. 


Landscape 

Mental 

Urban 

Elements 

C  oin  po  n  g  n  t 

Image 

Subsystem 

in  Art 

Strip 

Path 

Highway  strip 

Line 

Zone 

Edge 

Fringe  (urban-rural) 

Border 

Center 

Node 

Central  business  district 

Point  (foreground) 

Background 

District 

Neighborhood 

Background 
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the  example  of  natural  landscape  components,  shape  is  an  emergent  property 
of  the  interaction  of  ecosystem  and  topographic  processes. 

Comparison  of  Energy  in  Landforms  and  Ecosystems 

Embodied  energy,  actual  energy,  and  transformation  ratios  for  the  ten 
landforms  considered  are  given  in  Tables  33  and  34  on  an  area  and  volume 
basis,  respectively.     Transformation  ratios  are  the  same  or  very  nearly 
the  same  from  either  an  area  or  volume  view.     An  exception  is  the  glacier 
whose  great  thickness  (up  to  100  times  larger  than  the  other  landforms) 
causes  its  area-based  quality  factor  to  be  very  low,  though  the  volume- 
based  quality  factor  is  within  the  range  of  the  other  values.     The  glacier 
also  has  the  largest  actual  energy  due  to  its  large  mass  and  elevation, 
which  contributes  to  its  relatively  low  quality  factor. 

Landform  transformation  ratios  range  over  8  orders  of  magnitude  from 
spoil  mound  at  106  solar  Cal/Cal  to  coral  reef  at  lO1^  solar 
Cal/Cal.     These  values  are  higher  than  most  other  storages  (Odum  et  al . 
1978a)   that  have  been  calculated.     The  spoil  mound,   formed  by  mechanical 
energy  of  strip  mining,  has  lower  quality  than  landforms   formed  by  natural 
energies  over  long  development  times.     Of  the  natural  landforms,  the  two 
desert  examples,   sand  dune  and  arroyo,  have  the  smallest  transformation 
ratios  due  primarily  to  short  development  times  and  to  lesser  energy 
inputs.     Thus  the  spoil  mound  is  related  to  desert  conditions  as  will  be 
noted  concerning  vegetation  development.     Three  wetland  systems — palsa, 
salt  marsh,   and  seagrass-dominated  mud  mounds — have  similar  intermediate 
quality  factors,  on  the  order  of  109  to  1010  solar  Cal/Cal,   in  spite 
of  having  very  different  components  and  environmental  backgrounds. 


118 


Table  33.     Energy  transformation  ratios  of  land  forms  on  an  area  basis. 


Land  form 


Embodied 
Energy 
Input , 
109  solar  Cal/m2 


Actual 
Energy, 
Heat  Cal/m2 


Transformation 
Ratio, 

1010  solar  Cal/heat  Cal 


Floodplain 
Coral  reef 
Salt  marsh 
Glacier 
Palsa 
Sand  dune 
Arroyo 
Oyster  reef 
Mud  mound 
Spoil  mound 


1,000 
51,000,000 

10.4 
260,000 

9.2 
30 
5.1 
11,000 
75 
0.02 


1.6 

242.9 
5.  1 

272,000,000 

2.5 
70.9 
114.  7 
5.5 
2.4 
20.4 


62.5 
21,000 

0.  2 

0.0001 
0.4 
0.04 
0.004 
200 
3.1 

0.0001 
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Table  34.     Energy  tr 

ansformation  ratios 

of  land  forms 

on  a  volume  basis. 

Landform  10^ 

Embodied3 
Energy 
Input, 

cnl  a r    Ca 1  / 

Energy, 

Wo  ah            1  /m3 

tie  a  l  ua  i  /  m 

Trans  format  ion 
Ratio, 

iu      solar  Cal/heat  Cal 

Floodplain 

1,000 

1.6 

62.5 

Coral  reef 

7,300,000 

26 

28,000 

Salt  marsh 

5.2 

2.6 

0.2 

Glac  ier 

2,000 

2,053 

0.  1 

Palsa 

2.  3 

0.  5 

0.46 

Sand  dune 

5.4 

9.6 

0.05 

Arroyo 

0.6 

13.6 

0.004 

Oyster  reef 

5 , 500 

2. 1 

262 

Mud  mound 

75 

1.9 

3.9 

Spoil  mound 

0.007 

8.2 

0.00009 

aEmbodied  energy  input  divided  by  depth  or  thickness  of  landform,  except  for 
spoil  mound,  which  was  calculated  separately  in  the  footnote  to  Table  22. 


120 


Highest  quality  landforms  are  the  two  reef  systems,  producer  reef  (coral 
reef  with  algae)  at  10^  solar  Cal/Cal  and  consumer  reef  (oyster  reef) 
at  1012  along  with  the  floodplain  at  1012  solar  Cal/Cal.     All  of 
these  systems  have  complex  living  structures  of  reef  skeleton  or  tree 
roots  that  contribute  directly  to  the  development  of  the  landform.  The 
tropical  landforms  (coral  reef  and  seagrass  mounds)  have  larger  quality 
factors  than  the  cold  climate  landforms  (glacier  and  palsa)  but  these 
differences  are  not  large. 

For  additional  comparison,  energy  parameters  were  calculated  for  eco- 
systems of  the  case  study  area  where  possible.     Energy  in  ecosystems  was 
calculated  in  two  ways,   from  the  energy  in  biomass  multiplied  by  a  trans- 
formation ratio  (Table  35)   and,   for  those  systems  with  available  data, 
from  development  (succession)   time  multiplied  by  embodied  energy  input  as 
was  done  with  landforms  (Table  36).     Comparison  of  embodied  energies  of 
landform  and  ecosystems  may  indicate  which  has  more  control  in  the  land- 
scape. 

Comparing  biomass-based  embodied  energy  from  Table  35  and  landform 
values   from  Table  33,   all  natural  landforms  have  much  higher  embodied 
energy  than  corresponding  ecosystems,  ranging  from  107  to  10^  times 
greater.       The  spoil  mound  landform,  however,   is  one  order  of  magnitude 
less  than  the  biomass  value  for  embodied  energy  per  area.     This  is  prob- 
ably because  the  form  of  the  mound  is  built  by  machinery  with  less  energy 
embodied  than  most  natural  landforms  and  because  of  the  high  amount  of 
biomass  in  the  spoil  mound  forest. 

Comparisons  can  be  made  for  three  ecosystems  by  using  succession- 
based  embodied  energy  from  Table  36.  The  significant  feature  of  these 
values  is  that  they  are  greater  than  those  based  on  biomass  (Table  35). 
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Table  36.     Succession-based  embodied  energy  values  for  ecosystems. 


Ecosystem 

Succession 
Time ,  yr 

Embod  ied 
Energy 
Input , 
109  solar  Cal/m2-yr 

Embodied 
Energy, 
10y  solar  Cal/m2 

Floodplain 

480  a 

6,800 

3,200,000 

Sand  dune 

800b 

1.2 

960 

Spoil  mound 

50c 

0.07 

3.5 

aAverage  of  data  from  Fonda  (1974),   Shel ford  (1954),   and  Nixon  (1975). 
bChadwick  and  Dalke  1965. 
cThis  study. 
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This  suggests  that  energy  is  embodied  in  ecosystem  characteristics  other 
than  biomass,  such  as  diversity,  niche  structure,  or  other  emergent  prop- 
erties.   Biomass-based  ecosystem  values  depend  only  on  energy  in  biomass 
and  leave  off  energy  in  ecological  complexity.     The  difference  between  the 
succession-based  value  and  the  biomass-based  embodied  energy  gives  a  mea- 
sure for  the  energy  in  ecological  complexity  which  is  at  least  1000  times 
greater  than  biomass  on  an  area  basis. 

Comparing  succession-based  embodied  energy  per  area  with  landform 
values,  ecosystems  have  more  energy  than  landforms.     The  results  of  these 
comparisons  are  contradictory  with  those   found  with  biomass  values,   and  it 
is  not  clear  which  ecosystem  value  is  most  appropriate  for  comparison.  Do 
ecosystems  control  landforms?     It  is  well  known  that  vegetation  regulates 
erosion  and  deposition  of  sediments  and  soils.     Also  the  reef  skeletons 
are  directly  deposited  by  the  metabolism  of  the  reef  organisms.  However 
there  is  an  old  school  of  thought  in  plant  ecology  called  Physiographic 
Ecology  (Cowles  1901;  Braun  1916;  Gano  1917),  which  holds  that  landforms 
determine  ecosystems  together  with  other  factors  such  as  climate.     It  is 
suggested  by  the  energy  calculations  given  here  that  the  invisible  struc- 
ture of  the  ecosystem  rather  than  the  physical  biomass  determines  the 
landscape.     More  thinking  is  required  concerning  interpretation  of  these 
numbers.     For  example,  the  effect  of  turnover  times  of  biomass,   form,  and 
ecological  complexity  is  one  question  to  be  considered. 

Energy  Effect  of  Landforms  as  a 
Contribution  to  Human  Economies 


Like  air,  land  is  often  taken  for  granted  due  to  its  vast  extent  and 
integral   importance  to  human  activities.     However,   land  and  landform  is 


125 


the  platform  for  human  economies  and  natural  ecosystems.     As  a  platform, 
land  is  an  energy  input  to  the  processes  and  uses  that  take  place  upon  it. 
The  rate  of  flow  of  energy  from  landforms  is  the  rate  at  which  they  are 
being  used  up.     For  an  elevated  land  form  this  is  the  rate  of  erosion  and 
for  a  depression  it  is  the  rate  of  filling.     Before  human  disturbance, 
natural  vegetation  was  adapted  to  control  erosion  and  deposition,  as  part 
of  the  production  process  shown  in  Figure  1.     In  the  use  of  a  landscape, 
however,  humans  convert  natural  vegetation  to  other  uses  or  substitute 
less  mature  types  of  vegetation  such  as  agriculture  and  lawns.     They  also 
disturb  the  land  surface  itself  in  addition  to  vegetation.     These  actions 
tend  to  increase  erosion  and  deposition  and  thus  increase  the  rate  at 
which  landforms  are  used.     Humans  affect  the  material  balance  of  land- 
forms,  which  is  the  geomorphic  equivalent  of  mineral  cycling,  by  drawing 
down  storages  built  over  many  years  by  natural  processes.     This  observa- 
tion has  also  been  made  for  human  impact  on  geochemical  cycles  (Odum 
1972).     Examples  of  energy  flow  for  three  landforms  are  given  in  Table  37. 
The  undisturbed  ecosystem  uses  between  102  and  105  less  landform 
energy  than  human  land  uses. 

How  this  accelerated  energy  is  actually  used  by  human  land  uses  is 
not  entirely  clear.     The  energy  is  being  expended  but  whether  or  not  it  is 
embodied  in  the  land  use  is  an  open  question.     Some  kinds  of  landform 
attract  development  such  as  river  mouths  and   floodplains,  hilltops,  and 
others.     When  these  locations  actually  generate  more  economic  output  than 
adjacent  sites  on  different  landforms  then  the  energy  flow  from  use  of  the 
landform  probably  does  get  embodied  into  the  land  use  assets. 

As  a  start  at  assessing  the  national   impact  of  human  use  of  land- 
forms,  a  preliminary  calculation  was  made  including  the  loss  of  wetlands 
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Table  37.  Comparison  of  the  use  of  landform  energy  by  three  systems, 
assuming  a  unit  area  of  1  km2  and  unit  elevation  drop  of  1 
meter. 


System  Type 

Sediment  Yield, 
t/km2-yra 

Actual  Energy, 
Cal/yrb 

Mature  forest,  Hubbard  Brook 

3.3 

7.7 

Mountain  clear-cut,  West  Virginia 

302 

700 

Construction  site,  Baltimore, 

Maryland 

49,000 

114,000 

aBormann  and  Likens  1979. 

bActual  energy  =  (mass) (gravitat ional  acceleration) (height) . 
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and  the  filling  in  of  reservoirs.     These  calculations  are  summarized  in 
Table  38.     The  sum  of  these  two  losses  alone  is  equal   to  one-third  the 
rate  of  loss  of  fossil  fuels  from  U.S.  reserves  (Leibowitz  1979).  Other 
landforms  being  consumed  include  channelization  of  rivers  (loss  of  mean- 
ders),  leveling  of  topography  by  agriculture,  and  many  others.     To  what 
extent  these  losses  are  balanced  by  the  construction  of  farm  ponds,  dredge 
and  mining  spoil  mounds,  and  other  human  landforms  is  not  known,  but  it  is 
probably  less  due  to  increased  erosion  rates.     The  cost  of  restoring  the 
lost  landforms,  if  it  were  found  necessary  to  do  so,  would  undoubtedly  be 
high. 

Landscape  Energy  and  Aesthetics 

The  human  appreciation  of  an  object  for  its  beauty  is  its  aesthetic 
value.    Although  often  difficult  to  evaluate,  aesthetic  appeal  of  land- 
scapes is  well  known.     Attraction  of  visitors  to  parks  is  one  measurement 
of  this  value.     Calculations  in  this  dissertation  have  suggested  that  dif- 
ferent landforms  can  have  widely  different  embodied  energies.     They  also 
probably  have  different  aesthetic  value.     Is  the  beauty  of  landscape 
related  to  the  energy  embodied  in  it?     Some  of  the  best  known  parks  are 
those  with  very  large  physical  dimensions  that  are  visible  or  made  visible 
to  visitors  to  the  park  from  observation  points.     Examples  are  the  Grand 
Canyon,  Niagara  Falls,  and  ocean  beaches.     Aesthetic  value  may  be  propor- 
tional to  area  visible.     Human  eye  perhaps   integrates   the  energy  of  the 
large  area,  with  visual  appeal  related  to  the  total  amount  of  energy  that 
can  be  seen.     Many  techniques  have  been  used  to  measure  aesthetic  beauty 
in  environmental  evaluations  (Arthur  and  Boster  1976),  but   these  "tend  to 
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be  overly  complex,  in  addition  to  being  highly  subjective  and  value-laden" 
(Redding  1973).  Potential  exists  for  a  quantitative  theory  of  landscape 
aesthetics  based  on  correlation  of  embodied  energy  calculations  with  inde- 
pendent measures  of  human  appreciation  such  as  visitor  statistics  or  ques- 
tionaires.  Odum  (1962)  discusses  the  "ecological  roots  of  the  aesthetics" 
and  raises  the  hypothesis  that  selection  for  appreciating  landscapes  may 
be  linked  to  survival  value. 

Spoil  Mound  Succession 

Phosphate  spoil  mound  succession,   like  succession  throughout  the 
coastal  plain  on  most  fallow  agricultural  fields  (without  fire),  proceeds 
from  pioneer  weeds  and  old  field  species  through  pines  to  hardwoods.  With 
fire,  pines  persist  without  hardwoods  (Chapman  1932).     Phosphate  spoil 
succession  is  shown  diagramat ically  in  Figure  29.     Colonization  by  weedy 
forbs  begins  directly  after  spoil  mounds  are  deposited.     These  include 
Eupatorium  capillifolium,  Ambrosia  artemisifolia,  Phytolacca  rigida, 
Heterotheca  subatillaris ,  and  others.     At  this  early  stage,  vegetation  is 
very  sparse  and  the  landscape  is  somewhat  desertlike  in  appearance.  In 
fact,  vegetation  pictures  of  badlands  (Brown  1971)   look  similar  to  early 
stages  of  spoil  mound  succession.     This  was  also  noted  by  DeVall  and 
Devane  (1949)   in  an  early  description  of  phosphate  spoil  mound  conditions. 
Grasses  such  as  Rhynchelytrum  repens  and  Paspalum  notatum  increase  in 
importance  over  forbs  in  the  first  5—20  years  after  mining.     Depending  on 
seed  source  location,   trees  (Myrica  cerifera)   invade  after  5—10  years 
resulting  in  savannalike  formations.     Gradually  trees  coalesce,  forming  an 
open  forest  canopy  with  pines  dominating,   perhaps  with  scattered  hardwoods 


Figure  29.     Main  stages  in  phosphate  spoil  mound  succession.  Development 
time  is  about  50  years  to  mature  hardwood  forest. 
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like  wax  myrtle,  wild  cherry,  and  oak  species.    As  succession  proceeds 
control  shifts  to  hardwoods,  which  presumably  constitute  climatic  climax 
conditions . 

Details  of  changes  in  ecosystem  characteristics  found  in  this  study 
are  similar  to  literature  reports  on  other  ecosystems.     Perino  and  Risser 
(1972)  and  Ewel  (1977),  using  different  methods,  described  similar  general 
patterns  of  vegetation  cover  increase  through  succession.     Snedaker  (1970) 
and  Ewel  (1971)  described  solar  radiation  passage  through  vegetation  in 
tropical  succession  showing  a  similar  general   trend  as  reported  here  in 
terms  of  optical  density.     Ewel  (1971)  and  Byer  and  Weaver  (1977)  document 
increases  in  tree  height  over  successional  time. 

Several  studies  have  reported  aspects  of  litter  and  soil  organic 
matter  development  over  succession.     Tisdale  et  al .   (1966)   describe  suc- 
cessional increases  in  litter  depth,  and  Crocker  and  Major  (1955),  Ewel 
(1968),  and  Perino  and  Risser  (1972)  describe  litter  weight  increases. 
Likewise  Billings  (1938),  Smith  (1940),  and  Hanks  (1972)  demonstrate 
increases   in  soil  organic  matter  over  succession.     Billings   (1938)  also 
provides  data  on  soil  profile  development  over  time  in  succession. 

Maturity  of  Phosphate  Forests 

The  forest  development  on  50-year-old  sites  is  perhaps  the  most 
striking  aspect  of  spoil  mound  succession.     The  existence  of  forest  eco- 
systems on  old  spoil  mounds  demonstrates  the  potential  of  natural  succes- 
sion as  a  reclamation  alternative.     Determining  the  maturity  of  these 
forests  is  a  prime  goal  of  this  study. 
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Maturity  is  a  common-sense  concept  used  in  discussions  of  age, 
development,  and  organization  in  a  variety  of  contexts.  Quantifying 
maturity  of  systems  is  the  central  problem  in  succession  studies.  Mar- 
galef  (1963,  1968,  1974)  suggested  the  ratio  of  biomass  to  productivity  as 
an  index  of  maturity.     E.   P.   Odum  (1969)  provided  general  criteria  on  eco- 
system maturity  including  aspects  of  community  energetics,  community 
structure,   life  history,  nutrient  cycling  and  stability.     Whittaker  (1966) 
used  the  ratio  of  aboveground  biomass  to  aboveground  net  productivity 
(after  Margalef)   and  suggested  ratios  about  40-50  for  the  most  mature  for- 
ests in  the  Great  Smoky  Mountains  (see  also  Skeen  [1974]   for  an  applica- 
tion) .     Data  appropriate  for  comparison  with  these  indices  and  criteria 
are  lacking  in  this  study.     Instead,  comparisons  of  species  compositions 
based  on  literature  of  Florida  succession  along  with  comparisons  of  eco- 
system characteristics  between  the  50-year-old  mine  sites  and  an  undis- 
turbed hardwood  hammock  (San  Felasco)   are  used  for  assessing  maturity. 

The  oldest  spoil  mound  forests  studied  here  are  dominated  by  live 
oak,  Quercus  virginiana ,  with  several  other  hardwoods  (water  oak,  sweet- 
gum,  and  wild  cherry)  and  a  few  pines  (slash  pine).     At  the  Payne  Creek 
site  some  large  pine  stumps  are  found,   indicating  logging,  which  might 
explain  the  difference  in  basal  area  between  this  site  and  Christina. 
Water  oak  and  sweetgum  appear  to  be  reproducing,  with  saplings  in  the 
understory,  while  pine  and  live  oak  do  not.     These  forests  are  similar  in 
composition  to  the  live  oak  forests  described  by  Laessle  and  Monk  (1961), 
who  reported  lack  of  reproduction  by  live  oak  in  advanced  natural  forests. 
This  suggests  that  live  oak  is  declining  at  the  study  sites  and  transition 
to  another  system  is  taking  place.     Sellards  et  al .   (1915)  described 
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similar  forests  that  they  call  sandy  hammocks  and  Laessle  (1942)  described 
similar  stands  that  he  terms  xeric  hammock. 

Detailed  studies  of  Florida  succession  in  the  absence  of  fire  gener- 
ally show  a  progression  from  pines  to  xeric  hammocks  to  climax  mesic  ham- 
mocks (Laessle  1942;  Ansley  1952;  Quarterman  and  Keever  1962;  Monk  1968; 
Alexander  and  Crook  1974;  and  Veno  1975).     The  spoil  mound  forests  seem  to 
fit  in  the  transition  between  xeric  and  mesic  hammock  types  of  this 
scheme.     They  are  definately  post-pine  as  indicated  by  the  dominance  of 
hardwoods.     In  terms  of  species  composition  the  forests  might  best  be 
described  as  xeric  live  oak  hammocks.     The  apparent  decline  of  live  oak 
and  the  existence  and  reproduction  of  more  mesophytic  hardwoods  (sweetgura 
and  water  oak)   also  suggests  that  transition  is  occurring  to  mesic 
hammock.     This  is  supported  by  Monk's  (1967)  classification  of  the  species 
at  the  50-year-old  sites  as  climax  pioneers  that  "invade  successional  com- 
munities and  are  influential   in  the  conversion  to  climax." 

As  a  summary  of  all  succession  data,  Table  39  compares  the  50-year- 
old  phosphate  forest  with  San  Felasco.     The  percent  difference  between  the 
forests  is  shown  in  the  final  column.     The  San  Felasco  forest  has  higher 
levels  for  all  parameters  except  percent  soil  organic  matter.     The  phos- 
phate forest  has  23%  more  organic  matter  in  surface  soils  compared  with 
San  Felasco.    A  possible  explanation  may  be  greater  soil  respiration  at 
the  more  mature  San  Felasco  forest  drawing  down  on  the  storage  of  organic 
matter  in  the  soil.     Difference  for  all  other  ecosystem  parameters  range 
from  a  low  of  13%  for  transmittane  of  light  through  the  canopy  to  a  high 
of  69%  for  number  of  tree  species.     Most  differences  are  between  20%  and 
50%  greater  development  at  San  Felasco  compared  with  the  phosphate  for- 
est . 
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Table  39.     Comparison  of  succession  parameters  at  the  phosphate  forests 
with  an  undisturbed  climax  forest.     Percent  difference  is 
defined  as  the  value  of  the  phosphate  forest  divided  by  the 
value  of  San  Felasco  Hammock. 


Phosphate  San  Felasco 

Index  Forest  Hammock  %  Difference 


Basal  area  (ra2/ha)  28.5  35.4  19 

Density  (trees/ha)  368  602  39 

Optical  density  0.86  1.02  16 

Tree  height  (m)  19.2  25.8  26 

Litter  depth  (cm)  2.5  4.7  47 

Litter  weight   (g/m2)  1250.4  1594.2  22 

Percent  organic  matter 

in  soil  8.3  6.4  +  23 

Tree  species  number  5  16  69 


A  principal  difference  between  the  two  types  of  forest  is  the  diver- 
sity of  tree  species.     San  Felasco  has  more  species  and  also  a  more  even 
species  abundance  distribution.     These  distributions  are  shown  in  Figure 
30  in  terms  of  each  species  contribution  to  the  total  basal  area.  Signif- 
icant differences  between  the  two  forests  are  more  dominance  in  the  phos- 
phate forests  by  the  most  common  tree  species  (more  than  twice  as  much 
dominance  by  live  oak  compared  to  magnolia  or  sweetgura)   and  the  number  of 
rare  species  (seven  species  at  San  Felasco  with  less  than  1%  of  the  total 
basal  area  compared  to  the  phosphate  forest  with  no  species  contributing 
less  than  4%) . 

Another  difference  between  the  phosphate  forest  and  San  Felasco  ham- 
mock is  the  relationship  between  basal  area  and  density.     Although  the 
basal  area  of  the  two  forests  is  somewhat  close,  the  tree  density  at  San 
Felasco  Hammock  is  almost  twice  as  large  as  density  on  the  old  spoil 
mounds.    Moreover  the  average  tree  diameter  of  the  phosphate  forest  is 
greater  than  the  average  at  San  Felasco  (27.5  cm  versus  23.5  cm,  respec- 
tively).    In  other  words,  there  are  fewer  trees  in  the  spoil  mound  forests 
but  their  average  size  is  larger  compared  with  San  Felasco  Hammock.  To 
examine  this  observation  in  more  detail,  diameter  size  class  distributions 
of  the  two  forest  types  were  constructed.     These  were  drawn  as  age  class 
pyramids  since  size  is  often  used  as  a  surrogate  for  age  in  trees.  This 
is  another  view  of  hierarchy,   in  this  case  in  terms  of  age.     Figure  31 
shows  percent  basal  area  versus  size  class  (31a)   and  percent  of  the  number 
of  individuals  versus  size  class   (31b).     The  size  pyramid   for  basal  area 
of  San  Felasco  Hammock  is  almost  symmetrical  about  the  middle  size  class. 
For  the  phosphate   forest  a  somewhat  bimodal  distribution  was   found  with 
less  basal  area  in  the  intermediate  classes.     This  could  indicate  that 
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Figure  30.     Comparison  of  basal  area  distributions  by  species  rank  at  (a) 
the  50-year-old  spoil  mounds  and  (b)  at  San  Felasco. 
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selective  disturbance  occurred  in  the  past,  such  as  logging.     The  distri- 
bution may  be  smoothed  out  as  the  size  classes  progress  through  time.  The 
shapes  of  the  size  class  pyramids  of  numbers  are  similar  for  both  forest 
types.     However,  numbers  in  the  larger  size  classes  fall  off  more  rapidly 
for  San  Felasco  than  for  the  phosphate  forest.     San  Felasco  Hammock  has 
more  trees  in  the  small  classes  and   fewer  in  the  large  size  classes  com- 
pared with  the  phosphate  forest. 

Table  40  compares  the  forests  in  terms  of  community  structure  or  com- 
plexity indices.     Again,  the  percent  difference  is  shown  in  the  final  col- 
umn.    Structure  is  greater  at  San  Felasco  for  all  indices  in  this  case. 
The  Holdrige  complexity  index  gave  the  greatest  difference  between  the 
forests  at  61%.     The  three  diversity  indices  give  differences  of  14%  for 
the  Simpson  index  to  67%  for  species  per  1000.     The  prediction  of  21  spe- 
cies per  1000  at  San  Felasco  is  very  close  to  the  value  of  22  species 
found  by  Monk  (1960)   in  more  extensive  sampling  of  this   forest.  Monk 
(1967)  also  calculates  average  Shannon  diversities  for  162  stands  in 
north-central  Florida.     The  highest  average  is  for  mesic,  calcareous  cli- 
max southern  mixed  hardwoods  at  3.07  bits  per  individual,  with  a  maximum 
of  over  3.60.     The  value  reported  in  this  study  for  San  Felasco  is  higher 
than  the  average  for  mesic  hammocks  but  lower  than  the  maximum. 

Floristic  similarity  of  the   forests  was  calculated  with  Sorenson' s 
index  at  19%  with  two  species  in  common.     This  yields  a  floristic  differ- 
ence of  81%.     The  similarity  index  was  also  calculated  using  tree  species 
lists  from  other  Florida  hardwood  hammocks.     Included  were  south  Florida 
tropical  hammock  (Harper  1921),   central  Florida  sandy  hammock  (Sellards  et 
al.   1915),  southern  mixed  hardwoods  (regional  description;  Quarterman  and 
Keever  1962),  north-central  Florida  mixed  hardwoods  (Monk  1965),  coastal 
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Table  40.     Comparison  of  community  structure  indices  at  the  phosphate 

forests  and  San  Felasco  Hammock.     Percent  difference  is  defined 
as  the  value  of  the  phosphate  forest  divided  by  the  value  of 
San  Felasco  Hammock. 


Phosphate  San  Felasco 

Index  Forest  Hammock  %  Different 


Holdridge  Complexity 

Index  44.4  115.1  61 

Simpson  Diversity 

Index  0.75  0.87  14 

Shannon  Diversity 

Index  2.16  3.27  34 

Species  per  1000  7  21  67 

Sorenson  Similarity 

Index  —  —  81 
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live  oak  hammocks  (Laessle  and  Monk  1961),  and  north-central  Florida  mesic 
hammock  (Veno  1976).     The  species  lists  are  included  as  Appendix  VI  with 
the  partial  matrix  of  similarity  values  given  in  Table  41.     Most  values 
are  low,  especially  for  the  phosphate  forest.     This  is  due  to  the  paucity 
of  species  found  in  the  study  forests,  presumably  since  sample  areas  were 
small.     The  phosphate  forest  was  found  most  similar  to  the  central  Florida 
sandy  hammock  and  San  Felasco,  as  expected,  was  most  similar  to  mixed 
hardwoods  of  north-central  Florida. 

Based  on  comparison  with  San  Felasco  and  literature  reports,  although 
well  developed,  the  phosphate  forest  is  definately  immature  and  some  suc- 
cessional  changes  toward  climax  can  be  expected.     The  length  of  time 
required  to  reach  climax  stage  is  not  known.     Table  42  provides  data  on 
successional  time  for  other  forests  in  the  southeastern  United  States  for 
perspective.     The  range  is   from  50  to  300  years.     Succession  on  phosphate 
spoil  probably  falls  in  the  lower  end  of  this  range.     A  complicating  fac- 
tor, however,   is  the  composition  and  quantity  of  seeding  input  to  the 
phosphate  forests.    When  the  oldest  sites  were  first  mined  50  years  ago, 
the  surrounding  region  was  much  less  developed  than  at   present  and  there 
were  probably  many  areas  of  little  disturbed  vegetation  to  act  as  seed 
source  to  the  mined  areas.     Under  these  conditions  succession  proceeded 
quickly  as  indicated  by  parameters  such  as  the  large  basal  area  of  phos- 
phate forests.     Over  time,  however,   the  surrounding  vegetation  has  been 
replaced  with  other  land  uses  or  at  least  has  been  significantly  changed. 
Few,   if  any,   seed  sources   for  climax  species  exist  any  longer,  which  may 
cause  succession  on  mined  land  to  be  arrested  to  subclimax  conditions. 
Some  ecosystem  parameters  such  as  cover,    litter  weight,   and  basal  area  may 
reach  climax  levels.     Other  parameters,  such  as  species  diversity  and 
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Table  41.     Partial  similarity  matrix  comparing  species  composition  of  the 
phosphate  forest  and  San  Felasco  with  other  Florida  hardwood 
forests  (see  text  for  forest  references). 


Phosphate  San  Felasco 

Forest  Hammock 


Phosphate  forest 

1 

0.  19 

San  Felasco  Hammock 

2 

0.19 

N.C.   Florida  mesic  hammock 

3 

0.24 

0.21 

Coastal  live  oak  hammock 

4 

0.13 

0.31 

N.C.   Florida  mixed  hardwoods 

5 

0.09 

0.63 

Southern  mixed  hardwoods 

6 

0.17 

0.57 

Central  Florida  sandy  hammock 

7 

0.32 

0.60 

S.   Florida  tropical  hammock 

8 

0.10 

0.07 
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aspects  of  size  class  distributions  that  represent  top  levels  of  popula- 
tion hierarchy  may  never  reach  climax  development. 

Energy  Sources  and  Reclamation 

The  proportion  of  the  total  budget  due  to  reclamation  found  in  this 
study  (0.3%,  see  Figure  27,  page  105)   is  very  similar  to  the  estimate  of 
0.2%  for  coal  mining  given  by  Ballentine   (1976).     Although  reclamation  is 
not  a  large  portion  of  the  total  budget  of  the  phosphate  industry  accord- 
ing to  calculations  given  here,   there  is  concern  and  need  to  make  the  best 
use  of  energy  in  reclamation.  The  mining  system  can  be  viewed  as  an  order- 
disorder  cycle  as  shown  in  Figure  32.     Extraction  and  processing  of  phos- 
phate pulses  ordered  land  into  disordered  wastes.     Reclamation  attempts  to 
restore  the  wastes  back  to  a  more  organized  state. 

A  major  problem  for  those  involved  in  reclamation  is  deciding  what  to 
restore.     Calculations  and  measurements  in  this  dissertation  showed  that 
spoil  mounds  may  not  be  wastes  but  actually  valuable  storages.  Mining 
topography  has  energy  embodied  in  it  that  is  wasted  if  leveled  to  premin- 
ing  grade  as  is  currently  required  by  law.     Not  only  is  landform  energy 
lost  in  leveling  spoil  mounds  but   fuels  to  drive  the  machinery  are  also 
wasted . 

If  spoil  mounds  must  be  leveled  to  flat  surfaces  this  may  be  done 
more  energy  effectively  by  blowing  them  up  with  dynamite  rather  than  with 
large  earth-moving  machinery.     Placing  a  small  charge  at  the  base  of  a 
mound  would  retain  some  topography  but  still  cut  down  the  angle  of  the 
slopes  and  disperse  a  large  amount  of  spoil.     The  cost  of  this  alternative 
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is  $2053/mound  vs  $5400/mound  if  leveled  by  bulldozers  (see  footnotes 
below)  . 

Laws  requiring  reclamation  criteria  need  to  be  analyzed  on  an  energy 
basis.     Those  laws  requiring  the  leveling  of  spoil  mounds  may  be  doing 
more  harm  to  the  overall  system  than  good.     More  encouragement  for  finding 
ways  to  use  and  restore  mining  wastes  through  low  energy  ecological  engi- 
neering rather  than  high  energy  mechanical   technology  are  needed. 

Recommendations  for  Using  Succession  as 
Reclamation  for  Phosphate  Spoil  Mounds 

The  problem  with  natural  succession  as  a  reclamation  alternative  is 
that  it  is  slow.     Techniques  for  speeding  up  succession  rates  would  make 
the  alternative  more  attractive.     Subsidizing  sensitive  aspects  of  the 
ecosystem  development  process  that  regulate  its  rate  may  speed  up  or 
"abbreviate"  succession  (H.   T.   Odum  1962).     As  noted  by  Odum  (1962:69), 
"where  succession  involves  accumulation  or  consumption  of  quantities 
between  starting  and  final  climax,   the  success ional  stages  may  be  abbrevi- 
ated by  artificial  addition  or  removal  of  that  quantity."     Such  activities 
fall  in  the  discipline  of  ecological  engineering  (H.  T.   Odum  1962,  1971; 
H.   T.  Odum  et  al .   1963b),  which  involves  the  use  of  small  energies,  rela- 
tive to  the  natural  sources,   supplied  by  man  to  produce  large  effects  in 
the  resulting  ecosystem  patterns  and  processes. 

Cost  of  leveling  spoil  mounds:   to  remove  rock  a  charge  of  0.5  lbs/yd3 
is  suggested  (Crimmins  et  al .   1972).     Since  spoil  material   is  less  dense 
than  rock,  assume  1/3  of  this  charge  would  be  required  to  blow  up  a  spoil 
mound.     Cost  of  dynamite  is  about  $35/25  kg  (Gregory  1977).     Cost  to  blow 
up  an  average  spoil  mound  (4000  yd3)   is  then:  Dynamite  cost  =  1/3(0.5 
lbs/ydJ)(4000  yd J/mound) ($35/25  kg)(2.2  kg/ lb)  =  $2 , 053 . 3/mound . 
Bulldozer  costs  is  $542.2/400  yd3  or  $1.35/yd3  (Godfrey  1976),  making 
the  cost  to  level  a  spoil  mound  =  $1 . 35/yd3 (4000  yd3)  =  $5 , 400 /mound . 
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A  simulation  model  of  succession  was  constructed  to  test  certain 
strategies  for  the  use  of  succession  in  reclamation  of  phosphate  spoil 
mounds.     The  model  developed  in  this  study  is  compared  with  other  models 
of  upland  succession  in  Table  43.     The  most  detailed  model  is  that  of 
Shugart  and  West  (1977)   for  Appalachian  Mountain  hardwoods  with  33  species 
simulated.     The  phosphate  forest  succession  model   is  the  most  aggregated 
of  these  models,  but  all  have  been  usefully  applied  to  practical  problems. 
The  model  was  designed  as  a  management  tool  to  assess  the  degree  to  which 
succession  may  be  able  to  be  accelerated  as  evaluated  by  time  to  climax. 
The  role  of  shading  and  germination  controls  on  the  rate  and  transitions 
of  succession  were  included.     Other  important  features  such  as  soil  mois- 
ture, seed  storages,  and  mortalities  and  aspects  of  mineral  cycling  could 
be  added  to  make  the  model  more  complete,  though  the  additional  complexity 
might  not  necessarily  add  to  the  utility  of  the  model.     Even  though  highly 
aggregated,  the  minimodel  of  succession  used  here  was  validated  with 
comparison  of  predicted  and  observed  litter  values. 

The  exercise  of  constructing  the  model  was  itself  useful  in  pointing 
out  deficiences  in  available  data  on  ecosystem  development,  such  as  rates 
of  seeding  in  natural  communities  and  control  functions  on  germination  of 
dominant  species.  Perhaps  research  funds  should  be  allocated  to  investi- 
gate these  factors  for  phosphate  vegetation. 

Two  kinds  of  simulations  were  studied,   seeding  and  addition  of  lit- 
ter.    Among  the  most  important  aspects  of  succession  are  the  input  of 
seeds  and  the  role  of  animals   in  initiating  and  maintaining  succession. 
Without  seeds  or  immigration  no  ecosystem  can  develop.     The  advantage  of  a 
nearby  seed  source  was  observed  several   times  in  this  study,  apparently 
allowing  succession  to  proceed  faster.     In  several  cases,  study  sites  that 
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seemed  advanced  for  their  age  had  nearby  seed  sources.     Model  results  sup- 
port this  observation  in  that  increased  seeding  did  accelerate  succession. 
Sell  (1974)  and  Gunderson  (1976)  describe  the  effects  of  seeding  deficien- 
cies in  other  ecosystems.     The  design  of  mine  sites  with  small  scattered 
areas  left  untouched  as  seed  sources  may  then  accelerate  succession  in 
surrouning  mined  land.     After  an  area  of  mined  land  develops  enough  to  act 
as  a  seed  source,  previously  untouched  surrounding  areas  could  then  be 
mined,   later  to  be  seeded  from  the  ecosystems  on  old  mine  land.  Whether 
this  plan  could  be  justified  on  an  overall  economic  basis  needs  to  be 
analyzed. 

The  accretion  of  a  litter  layer  across  the  soil  surface  of  an  ecosys- 
tem is  another  important  factor  in  succession.     As  Olson  (1963:328)  states 
concerning  succession,  "some  biological  development  may  be  specifically 
delayed  until  litter  and  soil  humus  have  approached  near  their  steady 
state  values."    Perino  and  Risser  (1972)  provide  an  example  of  litter 
deficits  in  terms  of  biomass  in  early  succession.     Litter  additions  were 
simulated  here  under  two  model  structures,  with  and  without  a  light 
threshold  control  on  the  germination  of  hardwoods.     With  light  control  no 
advantage  to  litter  addition  was  found,  but  without  this  limiting  action 
succession  was  accelerated  by  30%.     This  great  difference   in  predictions 
suggests  the  model   is  sensitive  to  germination  control  and  thus  requires 
further  attention.     Field  germination  tests  may  provide  useful  information 
concerning  this  question.     One  problem  in  the  use  of  litter  additions  is 
finding  an  inexpensive  analog  of  true   forest  litter.     The  latter  is  a  sen- 
sitive substance  that  both  insulates  and  aerates  the  soil,  has  vertical 
stratification  of  characteristics,   and  has  high  embodied  energy  (in  the 
standard  simulation  run,  Figure  25a,   litter  took  50  years  to  reach  steady 
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state).    Wood  chips,  hay,  and  other  kinds  of  mulch  may  be  poor  substitutes 
for  natural  litter.     Because  of  this,   seeding  strategies  may  be  a  more 
efficient  approach  to  accelerating  phosphate  succession. 

Succession  can  be  accelerated  even  faster  if  some  idea  of  the  succes- 
sional  sequence  is  known.     The  strategy  consists  of  shortening  the  length 
of  each  stage  in  succession  by  seeding  advance  stages  at  each  step.  In 
other  words,  suppose  a  transition  between  two  successional  stages  takes  10 
years  with  natural  seeding  rates.     It  may  be  possible  to  cut  the  transi- 
tion time  greatly  by  increasing  the  seeding  rate  of  the  advanced  stage 
after  allowing  the  previous  stage  a  minimal  development  of  say  5  years. 
This  does  not  mean  jump  a  stage,  but  rather  accelerate  it.     This  may  be 
necessary  since  each  successional  stage  "prepares"  or  conditions  the  site 
for  the  next  stage.     To  illustrate  this  principle,  an  example  of  how  it 
might  be  applied  to  phosphate  spoil  mound  restoration  is  given  below. 

Previous  work  has  identified  the  main  stages  of  succession  on  phos- 
phate spoil  mounds  (Figure  29,  page  130  )  along  with  transition  times.  The 
general  sequence  is  for  forbs  and  grasses  to  colonize  the  mounds  within  1 
year.     Scattered  trees   first  appear  after  10  years  and  within  20  years 
tree  canopies  close  forming  initial  forest  structure.     Near-climax  hard- 
wood hammock  conditions  can  be  reached  after  50  years.     The  accelerated 
succession  principle  uses  this  known  sequence  to  speed  succession.  A 
possible  restoration  schedule  might  be  as  follows: 

1.   Plant  weed  seeds  such  as  dogfennel  the  first   few  years  to 
develop  initial  plant  cover,  which  generates   litter  and  soil 
organic  matter. 
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2.  Plant  early  successional  shrub  and  tree  seeds  like  wax  myrtle, 
willows,  elderberry,  and  other  species  after  3  years.  These 
are  fast  growing  and  should  quickly  develop  a  shrub  woodland. 

3.  Finally,   plant  seeds  of  climax  species  after  about  10  years  to 
grow  up  in  the  shade  of  the  shrub  woodland.     Growth  will  be 
slow  and  some  clearing  or  herbiciding  of  shrubs  may  be  neces- 
sary to  help  climax  seedlings  emerge.     Near  climax  conditions 
could  be  reached  in  25  years  perhaps,  cutting  the  natural  suc- 
cession time  in  half. 

The  simulations  of  succession  with  added  seeds  support  the  idea  of  serai 
subsidy  suggested  above. 


APPENDIX  I 

KEY  TO  THE  SYMBOLS  IN  THE  ENERGY  CIRCUIT  LANGUAGE 


Energy  source  (forcing  function),  source 
of  external  cause. 


Heat  sink,  outflow  of  used  energy 


Energy  interaction,  one  type  of  energy 
amplifies  energy  of  a  different  quality 
(usually  a  multiplier) 


Economic  transaction  and  price  function 


Storage  (state  variables) 


On-off  switch  (digital  actions) 


Group  symbols  (1)  autocatalytic  self- 
maintenance  units,   (2)  production  units, 
and  (3)  general  purpose  box  for  miscel- 
laneous subsystems. 
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APPENDIX  II 
CASE  STUDY  DATA  TABLES 
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APPENDIX  III 

SUCCESSION  MODEL  COEFFICIENT  CALCULATIONS  (see  Figure  20) 


J  =  incident  sunlight  =  1,460,000  Cal/m2-yr; 
R  =  remaining  sunlight  =  365,000  Cal/m2-yr; 
Ql  =  hardwood  bioraass  =  15,700  g/m2  (Lugo  et  al .  1978); 
Q2  =  pine  biomass  =  9,000  g/m2  (Wiegert  and  Monk  1972); 

Q3  *  grass  biomass  =  500  g/ra2  (average  value   from  Odum  1960;  Golley  1965); 
Q4  ■  litter  mass  =  1184  g/m2  (model  was  scaled  with  this  value  as 
characteristic  of  later  stages  of  succession,  see  Table  27); 
K^Qj^R  =  sunlight  used  by  hardwoods  =  1,095,000  Cal/m2-yr; 

Kj  =  (1,095,000  Cal/m2-yr)/[(15,700  g/m2 ) (365 , 000  Cal/m2-yr)]  =  1.91  x  10"4; 
L^2R  =  sunlight  used  by  pines  =  1,095.000  Cal/m2-yr; 

LX  -  (1,095,000  Cal/m2-yr)/[(9,000  g/m2 ) (365 , 000  Cal/m2-yr)]   =  3.33  x  10"4; 
M1Q3R  =  sunlight  used  by  grasses  =  1.095,000  Cal/m2-yr; 
Mx  =  (1,095,000  Cal/m2-yr)/[(500  g/m2) (365 , 000  Cal/m2-yr)]   =  6  x  10"3; 
K2QiR  =  hardwood  gross  primary  production  calculated  with  data  from  Lugo  et  al , 
(1978); 

Gross  production  =  (17.5  g/m2 -day) (270  days/yr;  9-month  growing  season)  = 
4725  g/m2-yr; 

K2  =  (4725  g/m2-yr)/[(15,700  g/m2) (365, 000  Cal/m2-yr)]  =  8.25  x  10~7; 
L2Q2R  =  pine  gross  primary  production  calculated  with  an  average  from  Bayley 

et  al.   (1974)  and  Dugger  et  al .  (1975); 
Gross  production  =  (6.75  g  C/m2-day)(2  g  organic  raatter/g  C)(270  day/yr; 
9-month  growing  season)  =  3645  g/m  -yr; 

L2  =  (3645  g/m2-yr)/[(9000  g/m2) (365 , 000  Cal/m2-yr)]   =  1.11  x  10"6; 
M2Q3R  =  grass  gross  primary  production  calculated  from  data  from  Golley  (1965); 
Gross  production  =  794  g/m  -yr; 

M2  =  (794  g/m2-yr)/[(500  g/m2 ) (365 , 000  Cal/m2-yr)  =4.35  x  10"6; 
K3Q1  =  hardwood  respiration  calculated  by  difference  as  gross  production  - 

litter  fall,  3656  g/m2-yr; 
K3  =  (3656  g/m2-yr)/(15,700  g/m2)  =  0.233; 

L3Q2  =  pine  respiration  calculated  by  difference  as  gross  production  - 

litter  fall,  3119  g/m2-yr; 
L3  =  (3119  g/m2- yr)/( 9000  g/m2)  =  0.347; 

M3Q3  =  grass  respiration  calculated  by  difference  as  gross  production  - 

litter  fall,  487  g/m2-yr; 
M3  =  (487  g/m2-yr)/(500  g/m2)  =  0.974; 

N^4  =  decomposition  rate  =  total   litter  fall  as  steady  state  =  1069  g/m' 
N:  =  (1069  g/m2-yr)/(1184  g/m2)  =  0.903; 
K4Q4  =  hardwood  litter  fall  from  Lugo  et  al .   (1978)  =  1069  z/mZ-vr: 
K4  =  (1069  g/m2-yr)/(15,700  g/ra2>  "  °-068; 

L4Q2  =  pine  litter  fall  from  Wiegert  and  Monk  (1972)  =  526  g/m2-yr; 
L4  =  <526  g/m2-yr)/(9000  g/m2)  =  0.058; 

M4Q4  =  grass  litter  fall  average  from  Odum  (1960)  and  Golley  (1965)  =  307  g/m2-vr; 
M4  =  (307  g/m2-yr)/(500  g/m2)  =  0.614; 

51  =  hardwood  seed  germination,  assumed  =0.1  g/m2-yr; 

52  =  pine  seed  germination,  assumed  =0.1  g/m2-yr; 

53  =  grass  seed  germination,  assumed  1  g/m2-yr. 
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APPENDIX  III 
SUCCESSION  SIMULATION  PROGRAM 


1  CLEflR  <5QQ> 

3  INPUT  "WHAT  IS  THE  LENGTH  OF  TIME  TO  RUN";  TO 

5  G0SU3  2900 

6  INPUT  "OUTPUT  TO  CRT<0>  OR  PRINTER<1V';  F2 

7  INPUT  "WHAT  IS  THE  DESCRIPTION  OF  RUN";  Bt 

8  IF  F2=8  THEN  PLOT  12 

9  IF  F2=l  THEN  PLOT  4  :  PRINT  B$ 

10  J=l.  46E6 

12  Tl=5?5 

13  T2=250 

14  73=400000 

15  Nl=  9029 
20  Kl=l.  91E-4 
30  K2=3.  25E-7 
40  K3=.  233 

SO  K4=.  OSS 
6a  Ll=3.  33E-4 
62  L2=l.  11E-6 
64  L3=.  34? 
66  L4=  053 

72  Ml=6E-3 

74  M2=4.  35E-6 
76  M3=.  974 

73  M4=.  614 

80  si=y 

83  S2=0 
86  S3=l 
88  E1=0 
90  E2=0 

92  E3=0 

94  Q1=0 

95  02=0 

96  Q3=0 

97  Q4=fj 

93  INPUT  "WHAT  IS  DT";  DT 
100  T=T+DT 

HO  R1=J/<1+K1*Q1> 

115  R2=R1/<1+L1+Q2> 

120  R3=R2/<1+N1*Q3> 

125  P1=K2+Q1*R1 

130  P2=L2+Q2*R2 

135  P3=ri2+Q3*R3 

160  IF  Q4>T1  GOTO  llOO 

170  Q1=Q1+DT+<P1-<K3+Q1>--::K4+Q1)  >+Sl 

175  IF  Q4>T2  GOTO  1120 
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135  02=G2+DT*<P2-<L2*02>-<L4*Q2> >+S2 

187  IF  T>9  GOTO  6000 

193  02=Q2+DT*<P3-<H3*03>-<ri4+03>>+S3 

195  Q4=G4+DT*< CK4*Q1+L4*Q2+M4*Q3>-<N1*Q4>) 

200  IF  F2=0  THEN  20© 

210  0=01/200  :  F=l  :  GOSUB  5000 

220  0=02/200  :  F=2  :  GOSUB  5000 

220  0=03/50  :  F=3  :  GOSUB  5000 

240  0=04/109  :  F=4  :  GOSUB  5000 

270  F1=F1+DT   :   IF-Rb»±-THEN  PRINT  P$  :  Fl=8 

2S0  P*=A$  f\-->\ 

290  POKE  40333,10 

2O0  PLOT  29  :  PLOT  IS 

210  PLOT  2 

220  PLOT  T  :  PLOT  01/200 
220  PLOT  255 
240  PLOT  29  :  PLOT  28 
250  PLOT  2 

260  PLOT  T  :  PLOT  02/200 
270  PLOT  255 
2S0  PLOT  29  :  PLOT  17 
298  PLOT  2 

408  PLOT  T  :  PLOT  03/50 
410  PLOT  255 
420  PLOT  29  :  PLOT  21 
428  PLOT  2 

440  PLOT  T  :  PLOT  04/100 
450  PLOT  255 
683  IF  T<T0  THEN  168 
703  PRINT  "FINAL  VALUES" 

713  PRINT  "Ql=";  01,  "Q2=";  02,  "Q2=";  02,  "04=";  04 
399  IF  F2=l  THEN  PLOT  5 

983  END 

1130  IF  R2>T2  GOTO  178 
1181  Si=.  1 

1132  IF  21=1  GOTO  173 
1135  G1=Q1+E1 
1187  21=1 
1110  GOTO  178 
1128  S2=.  1 

1122  IF  22=1  GOTO  135 
1125  Q2=Q2+E2 
1120  22=1 
1125  GOTO  185 
2338  REM 

2313  FOR  1=1  TO  3  :  fi*=A*+"  :"   :  NEXT 

2028  FOR  1=1  TO  8   :  Rl$=fll*+"  !"   :  NEXT 

2030  F*="1234567898" 
2343  RETURN 

5883  REM    PLOTTING  ROUTINE  FOR  PRINTER 

5813  Q=INTOV> 

5820  IF  1X2  THEN  0=2 

5820  IF  Q>79  THEN  0=79 

5348  Z*=MID-T<F$,F,  1) 


5959  P*=LEFT$<P*,  Q-1>+Z*+RIGHT*<P& 89-Q> 

5069  RETURN 

6999  IF  fl=l  GGTO  199 

6995  Q3=Q3+E3 

6919  R=l 

6915  GOTO  199 
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APPENDIX  V 
CALCULATIONS  FOR  PHOSPHATE  MINING  AND 
FERTILIZER  OPERATION  MODEL  (see  Fig.  27) 


(a)  Embodied  energy  costs  of  purchased  inputs. 

Mine  inputs  (Greene  1978):   Electricity  =  112  x  106  KWH/yr;  embodied  energy  = 
(112  x  106  KWH/yr)(860.5  Cal/KWH) (8000  solar  Cal/Cal)  =  7.7  x  1014  solar 
Cal/yr.     Labor  =  6.5  x  106  S/yr;  embodied  energy  =  (6.5  x  106  $/yr)(2.2  x 
107  solar  Cal/$)  =  1.4  x  10l4  solar  Cal/yr.     Fossil   fuel   (oil  and  gas)  = 
1.7  x  106  Gal/yr;  embodied  energy  =  (1.7  x  106  Gal/yr)(3.5  x  104  Cal/Gal) 
(3100  solar  Cal/Cal)  =  1.8  x  10T4d  solar  Cal/yr.     Capital  investment  =  8  x 
107  $/20  yr  =  4  x  106  $/yr;  embodied  energy  =  (4  x  106  $/yr)(2.2  x  107 
solar  Cal/$)  =  0.9  x  1014  solar  Cal/yr.     Total  mine  inputs  (x  1014  solar 
Cal/yr):   Electricity,    7.7;  Labor,    1.4;  Fossil  fuel,    1.8;  Capital  investment, 
0.9;  Total  =  11.8. 

Fertilizer  plant  inputs  (Phillips  1975):   Sulfuric  acid  =  2.75  tons  acid/ ton  fer- 
tilizer; embodied  energy  =  (2.75  tons  ac  id/ ton) (1.2  x  106  tons/yr)($35/ton 
acid)(2.2  x  107  solar  Cal/$)  =  2.5     1015  solar  Cal/yr.     Electricity  =  3.3 
KWH/ton;  embodied  energy  =  (3.3  KWH/ton)(1.2  x  106  tons/yr) (860 . 5  Cal/KWH) 
(8000  solar  Cal/Cal)  =  2.7  x  1013  solar  Cal/yr.     Labor  =  0.51  man-hours/ ton ; 
embodied  energy  =  (0.51  man-hours/ ton) (1 . 2  x  106  tons/yr) ($5/man-hour) (2 . 2  x 
107  solar  Cal/$)  =  6.7  x  1013  solar  Cal/yr.    Maintenance  =  $11.47/ton; 
embodied  energy  =  ($11 . 47/ ton) ( 1 . 2  x  106  tons/yr)(2.2  x  107  solar  Cal/$)  = 
3.0  x  1014  solar  Cal/yr.     Taxes,  interest,  overhead  =  $13.48/ton;  embodied 
energy  =  ($13. 48/ton) (1 . 2  x  106  tons/yr)(2.2  x  107  solar  Cal/$)  =  3.6  x 
1014  solar  Cal/yr.     Capital  investment  =  40.87  x  106  $/15  yr  =  2.73  x 
106  S/yr;  embodied  energy  =  (2.73  x  106  $/yr)(2.2  x  107  solr  Cal/$)  =  6.0 
x  1013  solar  Cal/yr.     Total  fertilizer  plant  inputs  (x  lO1^  solar 
Cal/yr):  Sulfuric  acid,   2.50;  Electricity,   0.03;  Labor,  0.07;  Maintenance,  0.30 
Taxes,  interest,  overhead,  0.40;   Capital   investment,  0.06;  Total  =  3.36. 

Total  industry  purchased  inputs  =  Mine  input  +  Fertilizer  plant   input  =  1.2  x 

10l3  solar  Cal/yr  +  3.4  x  lO1^  solar  Cal/yr  =  4.6  x  1015  solar 

Cal/yr. 

(b)  Embodied  energy  in  phosphate   fertilizer  sold  to  market.     Phosphoric  acid 
fertilizer  is  19%  phosphoric  acid  by  weight,   the  remaining  portions  being 
filler.     The  molecular  weight  of  phosphoric  acid,  H3PO4,   is  98  g/mole  of 
which  32%  is  phosphorus.     The  weight  of  phosphorus   in  fertilizer  is  (1  ton) 
(106  g/ton)(0. 19X0.32)  =  6.08  x  104  g  P/ton  fertilizer.     Using  the  embodied 
energy  per  gram  of  phosphorus,  7.63  x  10    solar  Cal/g,   from  Odum  and  Odum 
(1980)   and  the  annual   fertilizer  production,   the  total  embodied   flow  to 
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market  is  (1.2  x  106  tons/yr) (6. 1  x  1 04  g  P/ton)(7.63  x  105  solar  Cal/g)  = 
5.58  x  1016  solar  Cal/yr. 

(c)  Natural  energy  input.    Embodied  energy  f low  =  (2.1  x  10^  solar  Cal/m2,yr; 
see  Table  22)(400  acre/mine) (4047  m2/acre)  =  3.4  x  10l2  solar  Cal/yr. 

(d)  Embodied  energy  input  to  reclamation 

Clay  ponds.     2500  acres  of  clay  ponds  are  produced  each  year  in  the  central 
Florida  phosphate  district  (Anonymous  1974).     Dividing  by  the  average  annual 
production  of  34,000,000  tons  (Phosphate  Land  Reclamation  Study  Commission  1978) 
yields  (2500  acres/yr) /(34, 000, 000  tons/yr)  =  0.0001  acres/ton  rock.  Dollar 
cost  of  clay  pond  reclamation  is  $l,313/acre  (Phosphate  Land  Reclamation  Study 
Commission  1978),  making  the  embodied  energy  cost  per  ton  (0.0001  acres/ ton) 
($1313/acre)(2.2  x  107  solar  Cal/$)  =  2.9  x  106  solar  Cal/ton,  or  (2.9  x 
106  solar  Cal/ton  rock) (3. 31  ton  rock/ ton  fertilizer)  =  9.6  x  106  solar 
Cal/ton  fertilizer. 

Overburden  spoil.     An  average  mine  operates  on  400  acre/yr  (EPA  1978)  and 
extracts  5  x  10°  tons  rock/yr  (Greene  1978).     Acres  to  be  reclaimed  per  ton  of 
product,  assuming  onequarter  of  the  area  in  spoil  mounds,   is  (0.25)(400 
acres/yr) (5  x  10"  tons/yr)  =  0.00002  acres/ton.     Dollar  costs  of  spoil  recla- 
mation range  from  a  low  of  $300/acre  (Wang  et  al.   1974)   to  $1587.5/acre  (Phos- 
phate Land  Reclamation  Study  Commission  178).     An  intermediate  value  of 
$750/acre  is  used  here  as  suggested  by  Wange  et  al .   (1974).     Embodied  energy 
cost  of  reclamation  per  ton  is  (0.00002  acres/ ton) ($750/acre) (2 . 2  x  106  solar 
Cal/$)  =  3.3  x  105  solar  Cal/ton  rock,  or  (3.3  x  105  solar  Cal/ton  rock) 
(3.31  ton  rock/ ton  fertilizer)  =  1.1  x  10^  solar  Cal/ton  fertilizer. 

Sand  tailings.     Reclamation  costs  of  sand  tailings  are  given  in  terms  of  $/acre- 
ft.    The  amount  of  acre-ft  in  3.31  tons  of  sand  (total  sand  waste  per  ton  fer- 
tilizer, see  waste  volume  footnote  below)  =  [(1  ton  sand)(10^  g/ton)/(1.4 
g/cm3)l/[(2.832  x  104  cm3)/  (ft3)]  -  (7.14  x  105  cm3)/(2.832  x  104 
cm3/ft3)  =  25.2  ft3  in  terms  of  acreft  =  (25.2  ft3)/(43,560  ft3/acre- 
ft)  =  6  x  10~4  acre-ft/ton  fertilizer.     Dollar  cost  of  reclamation  is 
$63. 76/acre-ft  (Phosphate  Land  Reclamation  Study  Commission  1978).  Embodied 
energy  cost  of  sand  reclamation  is  (6  x  10~4  acre-ft / ton) ($63 . 76/acre-ft ) 
(2.2  x  107  solar  Cal/$)  =  8.4  x  105  solar  Cal/ton  fertilizer.     Total  recla- 
mation costs  (x  10°  solar  Cal/ton):  Clay  pond,  9.6;  Spoil,   1.1;  Sand  tailings, 
0.8;  Total  =  11.5.     Total  reclamation  cost  based  on  the  annual  production  of 
fertilizer  plants  is  (11.5  x  106  solar  Cal/ton)(1.2  x  106  tons/yr)  =  1.4  x 
1013  solar  Cal/yr. 


APPENDIX  VI 

SPECIES  LISTS  USED  IN  SIMILARITY  CALCULATIONS 


Southern  Mixed  Hardwoods 

(Quarter-man  and  Keever  1962) 

Fagus  grandif lora 
Quercus  alba 
Liquidambar  styracif lua 
Quercus  laurifolia 
Magnolia  grandif lora 
Quercus  nigra 
Carya  tomentosa 
Carya  glabra 
Pinus  taeda 
Quercus  f alcata 
Nyssa  sylvatica 
Ilex  opaca 
Pinus  glabra 
Acer  barbatum 
Quercus  michauxii 
Tilia  spp. 
Carya  pallida 
Fraxinus  americana 
Ostrya  virginiana 

Coastal  Live  Oak  Hammock 
(Laessle  and  Monk  1961) 

Quercus  virginiana 
Quercus  laurifolia 
Carya  glabra 
Ilex  opaca 

Osmanthus  americanus 
Sabal  palmetto 
Magnolia  grandif lora 
Persea  borbonia 
I  lex  vomitor ia 
Juniperus  silicicola 

South  Florida  Tropical  Hammock 
(Harper  1921) 

Bursera  simaruba 
Sabal  palmetto 
Coccolbis  laurifolia 
Quercus  virginiana 
Sideroxylon  foetidissimum 
Eugenia  consuf a 
Ficus  spp. 
Persea  borbonia 
Ichthyomethia  piscipula 
Lysiloma  bahamensis 


North-Central  Florida  Mixed 
Hardwoods  (Monk  1965) 

Acer  barbatum 
Acer  rubrum 
Carpinus  carol iniana 
Carya  spp . 
Celtis  laevigata 
Cornus  f lor ida 
Fraxinus  spp . 
Ilex  opaca 

Juniperus  virginiana 
Liquidambar  styracif lua 
Magnolia  grandif lora 
Nyssa  sylvatica 
Osmanthus  americanus 
Ostrya  virginiana 
Persea  borbonia 
Pinus  glabra 
Pinus  taeda 

North-Central  Florida  Mesic 
Hammock  (Veno  1976) 

Quercus  virginiana 
Quercus  hemisphaerica 
Quercus  myrtifolia 
Quercus  chapmani  i 
Pinus  palustr is 
Magnolia  grandif lora 
Vacc inium  arboreum 
Osmanthus  americana 
Lyonia  ferruginea 
Carya  glabra 
Liquidambar  styracif lua 
Myr ica  cerifera 

Central  Florida  Sandy  Hammock 
(Sellards  et  al.  1915) 

Quercus  laurifolia 
Magnolia  grandif lora 
Persea  borbonia 
Quercus  geminata 
Carya  glabra 
Cornus  f lor ida 
Quercus  virginiana 
Ilex  opaca 

Liquidambar  styracif lua 
Quercus  nigra 
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Appendix  IV.  (continued.) 


CarYa  SP-  Pinus  clausa 

Ficus  breuifolia  Sabal  palmetto 

Swietenia  mahogoni  Ostra  virginiana 

Sapindus  saponaria  Tilia  sp. 
Roystonea  regia 
Krugiodendron  ferreum 
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